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located  at  the  same  depths  below  the  water  surface  and  if 
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to  the  vertical  plane  containing  both  the  source  and  the 
receiver.  Experimental  results  -{3^  from  a  model  tank  operated 
under  various  wind  conditions  have  verified  the  existence  of 
unequal  Doppler  sidebands  under  these  conditions.  However, 
strong  asymmetries  in  sidebands  were  also  observed  under 
conditions  under  which  the  previous  theory  would  have  pre¬ 
dicted  no  asymmetries  at  all.  For  example,  strong  asymmetries 
have  been  observed  in  the  crosswind  direction. 

This  research  is  to  explain  the  asymmetries  not 
predicted  in  previous  theory  but  discovered  under  laboratory 
operations,  and  to  predict  other  asymmetries  not  yet  found 
in  either  the  experimental  or  the  theoretical  phase.  __ 

The  general  approach  to  the  solution  of  this  problem 
is  via  the  Kirchhoff  integral,  together  with  a  choice  of 
certain  surface  wave  models.  In  this  way,  it  is  shown  that 
asymmetries  can  arise  not  only  from  unequal  depths  of  source 
and  receiver,  but  also  from  small  misalignments  in  the 
measuring  system.  Furthermore,  it  is  shown  that  the  usually 
ignored  skew  in  the  surface-slope  probability  density  function 
also  contributes  to  the  observed  results. 
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Chapter  I.  Frequency  spreading 


1.0  Introduction 

When  an  acoustic  signal  Is  scattered  from  a  water  surface, 
the  received  signal  In  general  depends  on  the  characteristics  of 
the  reflecting  surface.  If  the  water  surface  is  a  moving  rough 
surface,  and  If  the  acoustic  signal  Is  a  pure  sinusoid,  a  substantial 
amount  of  frequency  spreading  around  the  transmitted  frequency  can  be 
expected  in  the  received  signal  [1]  -  [18].  This  frequency  spreading 
behavior  is  the  result  of  the  amplitude  and  phase  modulations  of  the 
acoustic  signal  because  of  surface  roughness.  Gulin  [36]  has 
described  the  amplitude  and  phase  modulations  of  a  sinusoidal  surface 
from  a  physical  optics  point  of  view.  Parkins  [4]  generalized  the 
sinusoidal  surface  case  to  a  time  varying  ocean  situation  where  he 
found  the  scattered  acoustic  spectrum  was  related  to  the  surface 
spectrum.  Scharf  and  Swarts  [S]  and  Eggen  [11]  have  developed  a 
simple  model  of  high  frequency  scattering  which  gives  a  fairly 
general  qualitative  description  of  scattered  field  spectral 
characteristics.  They  found  the  bandwidth  of  the  frequency  spreading 
function  to  be  a  function  of  surface  roughness.  Clay  and  Medwin  [7] 
have  looked  at  the  problem  in  the  time  domain  and  derived  the  temporal 
correlations  of  the  scattered  sound  field.  McDonald  and  Tuteur  [8] 
have  considered  the  surface  scattering  function  which  contains  both 
time  delay  and  frequency  spreading.  All  of  the  above  papers  use 
the  physical  optics  method  (Kirchhoff  integral*) .  Harper  and  Lablanca 
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[1],  [44],  and  Kuperman  [2],  [46]  have  considered  the  scattering 
problem  from  a  perturbation  point  of  view.  These  theoretical 
treatments  of  the  scattering  problem  discovered  a  spectral  asymmetry 
that  had  previously  escaped  notice. 

There  has  also  been  a  number  of  experimental  studies  of  the 
scattering  phenomenon.  Typical  are  sea  experiments  of  Roderick 
and  Cron  [16],  Brown  and  Frisk  [17],  and  model  tank  experiments 
such  as  those  of  Zornig  [3]  or  Gazanhes  et  al  [13],  etc.  There 
appears  to  have  been  no  single  theory  that  explains  all  the 
details  of  the  frequency  spreading  behavior  that  have  been  discovered 
in  these  experiments.  For  instance,  the  perturbation  method  of  Harper 
and  Labianca,  and  of  Kuperman  showed  that  the  frequency  spreading  function 
may  be  asymmetric  for  upwind  (or  downwind)  transmission,  but  it  does 
not  explain  the  observed  asymmetry  in  the  crosswind  direction.  Further¬ 
more,  it  applies  only  in  small  surface  roughness  situations.  For 
large  surface  roughness  the  physical  optics  method  has  generally  been 
used;  however,  if  the  usual  simplifying  approximations  are  applied, 
this  method  gives  no  spectral  asymmetry  at  all. 

f 

In  this  study,  the  frequency  spreading  function,  including  both 
the  amplitude  and  bandwidth,  is  examined.  The  object  of  this 
investigation  is  to  develop  a  model  and  a  formula  from  which  the 
general  frequency  spreading  function  can  be  predicted,  given  the, 
scattering  geometry  and  surface  conditions. 

The  analytical  solution  to  this  problem  is  based  on  the 
Kirchhoff  integral  method.  The  mathematical  results  presented  are 
for  arbitrary  scattering  geometry  where  the  locations  of  source  and 
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receiver  are  not  fixed.  These  results  are  specialized  for  upwind 
(or  downwind)  and  crosswind  geometries,  and  detailed  properties  of 
the  function  are  discussed  mainly  for  these  special  cases.  We  will 
find  that  the  important  surface  parameters  for  frequency  spreading 
are  the  second-order  moment  of  surface  waveheight  and  the  second 
and  third-order  moments  of  surface  slopes.  Although  the  effect  of 
the  surface  slopes  has  frequently  been  neglected  in  previous  surface 
scattering  studies  (e.g.  see  [4]-[ll] , [13]) ,  our  analysis  shows  that 
surface  slope  statistics  may  play  an  important  role  in  determining 
the  shape  of  the  frequency  spreading  function. 

One  important  aspect  of  frequency  spreading  is  the  sideband 
asymmetry  which  was  discovered  in  both  theory  (perturbation  method) 
and  experiment.  We  will  show  that  similar  results  can  be  obtained 
by  use  of  the  Kirchhoff  integral  method  if  certain  usually  neglected 
higher-order  expansion  terms  of  this  integral  formula  are  retained. 

The  organization  of  this  study  is  as  follows:  Chapter  I 
describes  1)  the  frequency  spreading  function  and  its  general  behavior, 
2)  a  brief  review  of  the  sideband  asymmetries  found  by  the  perturbation 
method,  and  3)  a  historical  background  review  in  the  areas  of  underwater 
acoustical  surface  scatter.  Chapter  II  contains  a  review  of  the 
Kirchhoff  integral  and  the  necessary  boundary  conditions  in  order  to 
provide  a  solution  to  the  surface  scattering  problem.  Chapter  III 
considers  a  deterministic  surface  model,  from  which  the  solutions  to 
the  frequency  spreading  are  obtained.  Chapter  IV  contains  a  discussion 
of  the  scattering  of  acoustic  waves  from  a  random  water  surface.  Both 
large  and  small  roughness  cases  are  considered.  Chapter  V  gives  a 


comparison  between  the  theory  and  the  experimental  results  obtained 
from  a  model  tank.  It  Is  shown  that  misalignments  In  source-receiver 
geometry  can  cause  significant  shift  In  the  frequency  spreading 
function.  In  Chapter  VI,  the  study  of  frequency  spreading  behavior 
of  acoustic  surface  scattering  Is  considered  from  a  physical  point  of 
view.  Chapter  VII  offers  a  summary  and  conclusions.  Appendices  A  to 
J  provide  the  necessary  mathematical  background  and  computational 
details  for  the  solutions  presented  in  Chapters  II  to  V. 


1.1  Motivations 


Consider  a  monochromatic  acoustic  wave  of  frequency 
scattered  from  a  moving  rough  water  surface.  The  received  signal 
spectrum  will  generally  be  spread  in  frequency  around  the  angular 
frequency  If  the  surface  roughness  is  small  and  if  the  surface 

deformation  is  roughly  periodic,  the  spreading  consists  of  a  fairly 
distinct  sequence  of  sideband  frequencies.  This  is  the  result  of 
phase  modulation  of  the  transmitted  sinusoid  by  the  moving  surface, 
and  the  separation  between  the  spectral  lines  at  the  receiving  point 
is  roughly  equal  to  the  surface  frequency  [5].  If  the  surface  is 
rough  and  confused,  the  sidebands  merge  together  into  a  more  or 
less  continuous  spreading  of  frequencies  around  the  transmitted 
frequency  uiQ. 

In  several  recent  papers  [1],  [44]-[45],  E.Y.  Harper  and 
F.M.  Lablanca  showed  that  if  the  source  and  receiver  are  at 
different  depths  the  upper+  and  lower  sidebands  may  have  different 
spectral  amplitudes.  A  similar  result  was  obtained  more  recently 
by  Kuperman  [2],  [46].  The  amplitude  ratio  depends  on  the  direction 
of  the  surface  wave  motion  with  respect  to  the  source-receiver 
geometry.  Specifically,  if  the  receiver  is  at  a  larger  depth  than 
the  source,  and  If  the  wind  causes  the  surface  wave  to  travel  away 
from  the  source  and  toward  the  receiver,  then  the  upper  sideband 
amplitudes  are  larger. 


The  upper  (lower)  sideband  refers  to  band  of  frequency  higher 
(lower)  than  the  transmitted  frequency  Uq. 
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Figure  1.1  Crosswind  source-receiver  geometry 


The  existence  of  sideband  asymmetries  has  been  observed  in 
a  series  of  experimental  results  obtained  in  a  model  tank  under 
simulated  wind  blown  surface  wave  conditions  £3J,  [12 J.  These 
measurements  showed  the  predicted  asymmetries,  but  they  also 
showed  asymmetry  in  the  frequency  spreading  function  when  the 
source  and  receiver  were  at  the  same  depths,  and  the  source-receiver 
geometry  was  crosswind.  For  this  arrangement  the  Harper,  Labianca, 
and  Kuperman  theory  would  have  predicted  symmetrical  sideband 
structures.  It  is  clear,  therefore,  that  there  are  features  of  the 
spreading  mechanism  that  are  not  yet  included  in  these  theories. 

The  observation  of  sideband  asymmetry  in  the  model  tank 
experiment  can  be  explained  heuristlcally  as  follows.  Consider  a 
crosswind  source-receiver  configuration  in  Figure  1.1.  Rays 
reflected  from  upwind  facets  have  upward  doppler  shifts,  whereas 
those  reflected  from  downwind  facets  have  downward  doppler  shifts 
[21].  In  a  wind  driven  water  surface  the  slope  distribution  at  the 
windward  side  is  different  from  that  at  the  leeward  side  [55] ,  [56] . 

A  sketch  of  a  wind  blown  surface  wavs  is  shown  in  Figure  1.2j  note 
that  the  surface  slopes  on  the  windward  side  are  shallower  than  those 
on  the  leeward  side  [70].  The  area  of  up-wind  slopes  capable  of 
reflecting  rays  from  the  source  to  the  receiver  with  upper  doppler 
shifts  is  larger  than  the  area  of  the  downwind  slopes.  Hence,  one 
would  expect  the  upper  sidebands  to  be  larger  than  the  lower  sidebands. 
This  can  be  seen  in  the  experimental  results  of  J.G.  Zornig  shown  in 
Figure  1.3.  This  frequency  spectrum  was  taken  in  a  crosswind  condition 
with  the  source  and  receiver  at  equal  depths.  The  frequency  of  the 
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acoustic  wave  was  256  KHz.  As  shown,  frequency  components  higher 
than  the  transmitted  frequency  (i.e.  the  upper  sideband)  are 
larger  in  amplitude  than  the  corresponding  components  below  the 
transmitted  frequency  (i.e.  the  lower  sideband). 

We  see  that  in  order  to  investigate  this  phenomenon 
analytically,  it  is  necessary  that  the  analysis  takes  into  account 
the  details  of  surface  slopes.  Eckart  [31],  Gulin  [36],  Parkins  [4], 
McDonald  [38],  and  many  other  authors  generally  neglected  the  slope 
terms  in  the  Kirchhoff  integral  by  assuming  the  slope  variations  to 
be  small.  This  implies  that  amplitude  and  phase  modulations  on  the 
acoustic  signal  are  due  entirely  to  the  up  and  down  motion  of  the 
surface.  On  the  other  hand,  Tolstoy  and  Clay  [69]  included  slope 
terms  in  their  version  of  the  Kirchhoff  integral.  But  because  of 
the  assumption  of  a  directional  source,  these  slope  terms  were  com¬ 
bined  into  an  average  slope  term  at  the  surface  point  from  which 
most  of  the  acoustic  energy  was  reflected.  As  will  be  shown  in 
subsequent  chapters  of  this  work,  all  of  these  approximations  eliminate 
one  of  the  mechanism  that  causes  sideband  asymmetry. 
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The  experimental  technique  actually  involves  the  transmission  of 
narrow  pulses,  and  the  recording  of  the  pulse  responses.  After 
Fourier  transformation,  the  response  to  any  desired  exciting 
frequency  can,  however,  be  obtained. 
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1.2  Historical  Background 


The  problem  of  scatter  of  acoustic  signals  from  a  rough 
air-water  interface  has  been  extensively  studied  during  the  last 
few  decades.  An  excellent  survey  of  the  literature  to  1958  was 
compiled  by  Lysanov  [19].  A  more  recent  and  more  general  survey 
can  be  found  in  the  excellent  paper  of  Fortuin  [20],  which  covers 
work  up  to  1969.  The  following  review  leans  heavily  on  Fortuin' s 
paper,  but  extends  it  to  include  the  decade  1969  to  1979. 

The  scattering  of  sound  by  a  sinusoidal  water-air  interface 
was  first  studied  by  Lord  Rayleigh  [21].  This  was  the  first  attempt 
to  solve  the  wave  equation  in  combination  with  a  boundary  condition 
in  sound  scatter  problems.  Rayleigh's  solution  was  an  intuitive 
approach  which  assumed  that  the  scattered  sound  field  from  a 
pressure-release  sinusoidal  surface  could  be  described  by  a  discrete 
set  of  plane  waves  travelling  away  from  the  surface.  This  simple 
approach  received  criticisms  [22]-[30]  from  several  authors,  such 
as  Uretsky  [22],  Meecham  [28],  Heaps  [30],  who  questioned  the 
validity  of  the  boundary  conditions  assumed  by  Rayleigh.  All  of 
these  critics  concluded,  that  "the  Rayleigh  method  is  indeed  incorrect 
in  the  way  the  boundary  conditions  are  used"  [20, p.  1213],  however, 
"for  [relatively]  smooth  surfaces,  the  method  produces  results  that 
do  not  disagree  more  with  experimental  data  than  do  other,  more 
rigorous,  solutions."  Uretsky  [22]  developed  a  method  which 
improves  upon  that  of  Rayleigh.  It  is  based  on  a  more  rigorous 
mathematical  approach  with  matrix  computation  and  handling  techniques 
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to  evaluate  the  reflection  coefficients  for  the  various  orders  of 
sound  reflection  from  a  sinusoidal  surfaces.  His  method  was 
summarized  by  Barnard  et  al.  [27]  and  used  to  make  comparisons 
with  experimental  results  obtained  in  a  model  tank.  Satisfactory 
agreement  between  theory  and  experiment  was  obtained.  Both  the 
Rayleigh  and  Uretsky  theories  apply  to  sinusoidal  surfaces  only. 

In  1953,  Eckart  [31]  applied  an  integral  formulation  known 
as  the  Helmholtz  integral  [32]  to  the  solution  of  sound  scattering 
from  rough  boundaries.  He  assumed  two  boundary  conditions.  One 
is  that  the  water-air  interface  is  a  pressure-release  surface;  the 
other  is  that  the  first  derivatives  of  the  Incoming  acoustic  wave 
and  the  scattered  fields  are  equal  at  the  surface.  This  is  a 
frequently  used  boundary  condition  called  the  Kirchhoff  approximation. 

The  use  of  the  Kirchhoff  approximation  Implies  that  the  surface  is 
locally  flat,  l.e.  the  radius  of  the  surface  curvature  is  much 
larger  than  the  acoustic  wavelength,  and  the  surface  is  free  of  shadowing 
and  multiple  reflections.  Eckart  derived  a  surface  scattering  cross- 
section  coefficient  in  terms  of  the  incident  acoustic  frequency  and 
surface  spatial  spectrum.  The  Kirchhoff  integral  approach  (Helmholtz 
integral  plus  Kirchhoff  boundary  condition)  was  criticized  by 
Meecham  [33]  and  Mintzer  [34]  because  the  limit  of  the  validity  of  the 
Kirchhoff  boundary  condition  was  not  properly  considered.  It  has  in 
fact  been  said  [6]  that  Eckart  "had  obtained  significant  results  with 
minimum  mathematical  complexity  by  relying  on  a  highly  developed 
physical  insight  into  the  problem”.  Nonetheless  the  Kirchhoff  integral 
method  has  been  widely  used  in  studies  of  surface  scatter,  and  experi- 
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mental  verification  of  some  results  [36]  seems  to  indicate  that 
the  error  introduced  by  the  Kirchhoff  approximation  may  not  be  too 
serious.  The  formulation  used  by  Eckart  is  equivalent  to  the 
Fraunhofer  approximation  used  in  optics.  An  improvement  of  this 
approximation  is  the  Fresnel  or  second-order  approximation,  first 
applied  to  surface-scatter  problems  by  Feinstein  [35].  It  has 
also  been  used  by  Gulin  [36],  Melton  and  Horton  [37],  Medwin  and 
Clay  [7],  McDonald  and  Tuteur  [8],  etc.  The  Fresnel  corrected 
Kirchhoff  integral  is  usually  applied  where  the  illuminated  surface 
area  is  large.  It  has  many  applications,  such  as  the  study  of 
amplitude  and  phase  fluctuations  [12],  the  correlation  and  power 
distribution  of  scattered  sound  [7],  [38],  and  the  studies  of 
frequency  spreading  and  shift  in  forward  and  backward  scatters  [10]. 

A  different  approach  to  the  solution  of  surface  scatter 
problems  is  the  perturbation  method  introduced  by  Isakovich  [41] 
in  1957.  Its  original  application  was  to  solve  the  normal-mode 
transformation  problem  in  an  irregular  waveguide.  Bezrodnyi  and 
Fuks  [42]  used  it  in  the  solution  to  the  problem  of  amplitude  and 
phase  fluctuations  in  a  waveguide.  Wait  [43]  analyzed  the  reflection 
from  two-dimensional  periodic  sea  waves.  In  1975,  Harper  and 
Labianca  [44]  used  the  perturbation  method  in  the  study  of  the  spectral 
behavior  of  the  surface  scattered  acoustic  signal.  Their  studies  show 
the  existence  of  asymmetrical  sidebands  in  surface  scattering.  This 
result  was  corroborated  by  the  theoretical  study  of  Kuperman  [2] 
using  the  same  method,  and  by  experiment  [3].  The  perturbation  method 
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is  mathematically  more  rigorous  than  the  Kirchhoff  integral  method 
and  permits  fairly  precise  bounds  to  be  placed  on  the  validity  of 
the  simplifying  assumptions  [45].  However,  its  applications  are 
limited  to  problems  with  very  small  roughness. 

Both  the  Fresnel-corrected  Kirchhoff  integral  and  the 
perturbation  method  require  a  solution  to  the  wave  equation  subject 
to  boundary  conditions  at  the  water-air  interface.  If  the  form  of 
boundary  is  very  complicated,  it  can  be  extremely  difficult,  if  not 
impossible,  to  reach  such  a  solution.  Middleton  [47]  developed  a 
so-called  Quasiphenomenological  approach  which  eliminated  the  need 
to  simultaneously  satisfy  the  wave  equation  and  boundary  conditions. 

It  introduces  the  surface  (boundary)  in  terms  of  a  group  of 
randomly  distributed  point  scatterers.  Each  of  these  point  scatterers 
has  its  own  system  function  and  directivity  pattern.  This  gives  the 
model  the  ability  to  handle  problems  of  general  scatter  geometry  with 
very  complex  boundary  conditions.  However,  application  of  this  very 
general  theory  is  severely  restricted  in  practice  by  the  difficulty 
of  relating  measured  surface  parameters  to  rtatistical  behavior  of 
the  scatterers. 

The  commonly  used  surface  parameters  are  surface  waveheight 
and  surface  slopes.  In  most  of  the  analyses  that  require  descriptions 
of  the  water  surface  boundary,  it  is  assumed,  with  exceptons  e.g. 
[49],  [SO],  that  the  water  surface  height  and  slopes  are  stationary 
Gaussian  processes  [4]-[5],  etc.  Experiments  performed  at  sea  and 
in  model  tanks  have  shown  that  this,  in  fact,  is  close  to  the  real 
situation.  For  example,  Kinsman  [51]  used  a  capacitance  probe  held 
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at  the  water  surface  to  record  the  surface  displacement.  The 
result  of  his  measurement  indicated  the  surface  height  distribution 
was  close  to  Gaussian  with  zero  mean.'*'  This  observation  was 
corroborated  by  Weissman  [52]  in  his  ocean  waveheight  measurement 
with  a  two  frequency  radar  interferometer  and  by  Zomig  [53]  in  his 
height  measurement  in  a  model  tank.  However,  Spindel  [54]  found 
that  the  joint  distribution  of  two  points  on  the  surface  is 
occasionally  far  from  Gaussian.  In  particular,  he  found  that  the 
conditional  distribution  of  the  surface  at  one  point,  given  its 
value  at  the  second  point,  was  occasionally  multimodal,  especially 
at  high  winds.  Cox  and  Munk  [55]  studied  the  slope  distribution  of 
a  wind-driven  surface  using  aerial  photograph  technique.  Their 
method  consisted  "in  photographing  from  a  plane  the  sun's  glitter 
pattern  on  the  sea  surface,  and  translating  the  statistics  of  the 
glitter  into  the  statistics  of  the  slope  distribution".  Figure  1.4 
contains  the  result  of  the  slope  measurement  by  Cox  and  Munk.  The 
top  plot  was  the  slopes  measured  crosswind  and  had  a  symmetrical  structure 
with  respect  to  zero  slopes.  The  bottom  plot  was  the  up- down  wind 
slope  distribution  which  showed  skewness  to  the  upwind  side,  i.e.  the 
most  probable  slope  for  the  sea  surface  was  negative.  This  distortion 
probably  was  due  to  a  wind  stress  effect,  as  remarked  by  Kinsman  [51], 
p.  350].  Schooley  [56]  made  the  same  sea  slope  measurement  by  a 
similar  optic  method  except  the  sun  was  replaced  by  a  flashing  gun  attached 
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The  flat  surface  level  is  chosen  as  the  reference  zero. 


Figure  1.4  The  probability  distribution 
of  surface  slopes.  The  upper  curves  are 
for  crosswind;  the  lower  curves  are  for 
downwind.  The  solid  curves  refer  to  the 
observed  distribution,  the  dashed  to  a 
Gaussian  distribution  of  equal  mean 
square  slope  components.  The  thin 
vertical  lines  show  the  scale  for  the 
standardized  slope  components 

—  -ir'  for  crosswind  and 
o  9x 
c 


—  for  downwind. 

°D 

The  wind  direction  is  in  the  y  axis; 
oc  and  Op  are  the  r.m.s.  slope  in  x  and 
y  directions  respectively.  The  heavy 
vertical  segments  show  the  surface  slope 

angle  for  5°,  10°, . ,25°.  In  the 

lower  curve,  the  negative  slopes  are  on 
the  upwind  side  and  the  positive  slopes 
are  on  the  downwind  side. 


(From  Cox  and  Munk,  1954) 
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to  the  camera.  The  same  skewness  in  up-down  wind  slopes  has  been 
observed  by  Wu  [57]  and  Zornig  [58]  in  model  tank  experiments. 

There  are  two  reasons,  as  pointed  out  by  Kinsman,  why  the 
waveheight  and  slopes  are  not  strictly  Gaussian.  One  is  that  the 
wave  motion  must  satisfy  the  Bernoulli  equation  and  the  free-surface 
boundary  conditions.  The  other  is  that  waves  are  limited  in  height 
by  the  wave  breaking  phenomenon.  In  summary,  it  appears  that  even 
though  the  surface  distribution  is  generally  not  strictly  Gaussian, 
the  assumption  of  a  Gaussian  distribution  is  not  unreasonable,  and  it 
has  been  widely  used. 

Another  factor  to  be  considered  in  the  development  of  a 
scattering  theory  is  shadowing.  Shadowing  occurs  when  a  part  of  the 
surface  is  screened  by  some  other  parts  of  the  surface.  In  surface 
scatter  terms,  this  means  the  illumination  of  incoming  radiation  at 
these  parts  of  the  surface  is  interrupted.  Shadowing  is  more 
significant  when  the  acoustic  source  and  receiver  are  very  close  to 
the  scattering  surface.  Neglecting  this  shadowing  effect  may  cause 
large  errors  if  the  theory  is  based  on  the  assumption  of  total 
illumination  on  the  water  surface,  e.g.  the  Kirchhoff  integral 
method.  Beckman  [59]  proposed  a  simple  way  to  treat  the  shadowing 
by  introducing  a  shadowing  function  S  which  equals  unity  on  the 
illuminated  parts  of  the  surface  and  equals  zero  on  the  screened  parts. 
In  other  words,  the  shadowing  function  S  is  the  probability  that  this 
part  of  surface  will  get  illuminated.  He  applied  this  concept  of  a 
shadowing  function  to  the  study  of  backseat ter ing  from  a  composite 
rough  surfaces  [60].  However,  a  computer  simulated  experiment  on 
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shadowing  fn  backscatter  by  Brockelman  and  Hagfors  [61]  showed 
marked  disagreement  between  their  results  and  those  of  Beckmann. 

With  a  modification  to  the  Beckmann  theory  that  includes  the 
effect  of  slopes,  Wagner  [62]  developed  formulas  which  consider 
geometric  shadowing.  A  comparison  between  his  results  and  those  of 
experiment  showed  an  excellent  agreement  between  theory  and  experi¬ 
ment  [20].  Lynch  and  Wagner  [63],  [64]  used  this  geometric  shadowing 
to  correct  the  energy  loss  in  high  frequency  scatter  from  random 
rough  surfaces.  A  simplified  method  to  evaluate  Wagner's  shadowing 
function  was  given  by  Smith  [65].  He  omitted  the  correlation  between 
surface  height  and  slopes  but  obtained  a  result  which  was  not  very 
different  from  the  complete  solution  obtained  from  Wagner's  method. 
Gardner  [66]  used  Wagner's  theory  to  explain  near  grazing  backscatter 
through  the  Kirchhoff  integral  approach.  However,  it  was  pointed  out 
in  the  experiment  by  Novarini  and  Medwin  [67]  that  Wagner's  theory 
of  geometric  shadowing  is  inadequate  in  the  sea  experiment  because 
of  the  proven  existence  of  bubbles  below  the  sea  surface  [68].  In 
our  analysis  of  the  scattering  problem,  shadowing  has  been  ignored. 
The  analysis  therefore  is  limited  to  relatively  large  grazing  angles. 
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Chapter  II.  Mathematical  preliminaries 
2.0  Introduction 

In  this  chapter,  we  introduce  some  basic  mathematical  concepts 
which  are  of  interest  in  this  analysis »  These  include  the  derivation 
of  the  Kirchhoff  integral  formula  to  solve  the  surface  scatter  problem, 
the  Fresnel  approximation  and  other  expansions,  the  definition  of 
surface  roughness,  and  the  beam  pattern  function  for  both  source  and 
receiver.  Most  of  the  material  presented  in  this  chapter  can  be 
found  elsewhere.  For  example.  Baker  and  Copson  [32]  and  Eckart  [31] 
give  a  derivation  of  Kirchhoff  integrals;  the  Fresnel  expansion  is 
a  second  order  approximation  and  is  commonly  presented  in  optics 
text 8  [71];  the  discussion  of  the  beam  pattern  function  (directivity) 
for  an  acoustic  source  can  be  found  in  texts  of  acoustics,  such  as 
Beranek  [72]. 

We  shall  regard  the  underwater  acoustic  surface  scattering 
process  as  linear  and  slowly  time-varying.  This  permits  us  to  use 
a  monochromatic  signal  source. 
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2.1  The  Kirchhoff  Integral 


Consider  a  function  p(r,  ioQ,  t)  which  satisfies  the  wave 
equation  within  a  volume  V  bounded  by  a  surface  S.  For  our  purpose 
p  is  regarded  as  the  complex  time  varying  amplitude  of  the  pressure 
at  a  location  defined  by  and  at  the  frequency  Uq  and  time  t. 
Assume  that  p  and  its  first-order  derivatives  are  finite  and 
continuous  inside  the  volume  V.  Then,  from  the  Helmholtz  theorem 
(see  Appendix  A) ,  the  value  of  p  at  a  particular  point  P  inside  V 


and  at  time  t  is  given  by 
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where  pg  is  the  pressure  amplitude  at  the  surface  S,  r^  is  the 

scalar  distance  from  the  surface  element  ds  to  the  point  P  as  shown 

in  Figure  2.1,  and  n  is  the  surface  normal  defined  to  be  polsitve 

in  the  outward  direction.  The  surface  pressure  p  and  its  normal 
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derivative  are  evaluated  at  the  time  t  -  ,  where  c  is  the  sound 

c 

velocity;  however,  for  the  time  scale  of  Interest  in  our  work  this 
retardation  delay  can  be  ignored  [74]. 

In  much  of  our  work  we  are  interested  in  sound  transmitted 
and  received  by  directional  hydrophones  with  directional  beam 
patterns.  The  effect  of  such  a  beam  pattern  is  to  give  unequal 
weights  to  different  parts  of  the  surface  S,  and  this  effect  can  be 
included  in  equation  (2.1.1)  by  changing  the  Green's  function 


j  ~rl 


J  — rl 


to  B. 


This  results  In 


"0 

ri 


*> '  if  ff'\ 


3p  Jc  1 

if-  -  ».  - >  1  •  di  <2'1'la> 


where  the  beam  pattern  function  BR  is,  strictly  speaking,  evaluated 

rl 

at  the  retarded  time  t - . 
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We  now  assume  that  the  water  surface  is  illuminated  by  a 
monochromatic  acoustic  source  which  has  a  beam  pattern  function  Bg. 
The  incident  acoustic  pressure  p^  measured  on  the  water  surface  is 
then  given  by 
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(2.1.2) 


where  rQ  is  the  distance  from  the  source  point  to  the  surface  element 
ds,  and  Pq  is  the  amplitude  of  the  acoustic  source. 

The  pressure  amplitude  p  appearing  in  equation  (2.1.1)  is 

s 

regarded  as  the  surface  reflected  pressure.  It  is  related  to  the 
incident  pressure  through  a  surface  boundary  condition.  For  the 
air-water  interface  the  surface  can  to  a  very  good  approximation  be 
regarded  as  being  a  pressure-release  surface;  hence  the  boundary 
condition  is 


p„  +  pi 


(2.1.3) 


Also,  the  normal  derivative  can  be  related  to  the  normal 

derivative  of  the  Incident  pressure  by  invoking  the  Kirchhoff  boundar 
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condition: 


(2.1.4) 


The  Kirchhoff  condition  implies  that  the  surface  is  "locally 
flat",  i.e.  the  radius  of  the  surface  curvature  is  much  larger 
than  the  acoustic  wavelength.  It  also  implies  no  shadowing  or 
multiple  reflections  on  the  surface,  i.e.  the  surface  is  totally 
illuminated.  As  briefly  discussed  in  Chapter  I,  there  are  criticisms 
about  the  validity  of  the  Kirchhoff  condition.  A  summary  of  these 
criticisms  can  be  found  in  [20],  [74].  However,  by  employing  these 
two  boundary  conditions  and  replacing  with  its  equivalent  from 
equation  (2.1.2),  one  can  convert  the  integral  in  equation  (2.1.1a) 


to  the  form 


3  ir<Vri) 


■  b  fl%  kk"*  ~  tt —  }•  ^  (2-1-5) 

8  0  1 


If  we  define  the  channel  transfer  function  H(a>Q,t)  as  the 
ratio  of  p(P,t)  and  p^,  equation  (2.1.5)  can  be  rewritten  as 

u0 

H<“o’t>  ‘  b  JJi*  {"A  -  •  di  < 

s  —  0  1 


(2.1.6) 


Thus  H(uQ,t)  represents  the  instantaneous  amplitude  and  phase  at 
the  receiver  due  to  a  unit  amplitude  sinusoidal  signal  with  frequency 
Wq  at  the  source. 
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I 

2.2  The  Fresnel  approximation  and  other  expansion  formulae 

» 

A  commonly  used  coordinate  system  in  surface  scattering 

studies  is  shown  in  Figure  2.1.  The  X-Y  plane  is  the  plane  of 

^  the  smooth  (flat)  surface.  The  positive  direction  for  the  Z  axis 

is  downward  into  the  water.  The  origin  of  coordinates  is  taken  to 

be  at  the  smooth-surface  specular  point.  The  source,  receiver  and 

,  specular  point  lie  in  a  plane  which  makes  an  angle  with  the 

{  K 

-Y  axis.  Both  source  and  receiver  have  grazing  angles  i|>0  with  the 
flat  water  surface.  The  distances  from  the  origin  to  source  and 
,  receiver  are  respectively  Tqq  and  r^ .  The  locations  of  source  and 

receiver  are  described  as  (XT ,  YT>  ZT)  and  (Xp,  YR,  ZR) ,  respectively, 
where 

-r00  cos  *0  cos  *R  XR  -  r10  cos  cos  *R 

r0o  cos  *0  8in  *R  YR  •  "r10  COS  *0  COS  *R  (2’2’1) 

r00  8lQ  *0  r10  Sin  *0 

We  further  define  £(x,y,t)  as  the  surface  displacement  at  a  point 
(x,y)  and  time  t,  with  positive  5  being  upward.  Therefore;  the 
)  distances  from  the  source  and  the  receiver  to  an  arbitrary  point  on 

the  surface  can  be  represented  by  Tq  and  r^  respectively,  where 

rQ  -  [(X-Xt)2+  (Y-Yt)2+  fe+Zj)2]1*  (2.2.2a) 

rx  -  [(X-X^)2  +  (Y-Yr) 2  +  C+Z^2]’5  (2.2.2b) 

t 
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The  integration  of  equation  (2.1.6)  cannot  be  carried  out  in 

closed  form  if  we  substitute  equations  (2.2.2a)  and  (2.2,2b)  directly. 

We  therefore  first  rewrite  Tq  and  r^  in  the  form 

r0  “  r00[1 - V2XtX  +  2YTY  ~  2ZT  C  -X2-'*2-C2)]is  (2.2.3a) 

■  r00 

rl  =  r10[1 - +  2YRY  ~  2ZR5  (2.2.3b) 

rio 

If  the  distances  rQQ  and  r^Q  are  large  compared  to  the 
dimensions  of  the  insonified  area,  i.e.  Tqq,  r^Q  »  X,Y,S,  we  can 
then  expand  equations  (2.2.3a)  and  (2.2.3b)  in  a  power  series  by 
use  of  the  binomial  expansion  formula 

(l+a)n  -  1  +  na  +  a2  +  -n(n~1>iS.-2).  a3  + .  (2.2.4) 

for  a2  <  1. 

We  further  assume  that  the  surface  .deformation  is  small 
compared  with  other  dimensions,  i.e. 

4(X,Y)  «  X,Y  «  r00,  r1Q  (2.2.5) 

Therefore  we  retain  the  first  order  term  in  4  and  the  X,Y  terms  up 
to  the  second  order  in  the  expansion.  The  result  is 


equation  (2.2.7)  can  be  simplified  with  the  aid  of  equation  (2.2.1) 
to  give 


+  r. 


00 


+  r10  + 


2c  sin  ij> 


O’ 


(2.2.8) 


Under  these  conditions,  the  path  dependence  on  X  and  Y  is  completely 
eliminated,  and  the  total  path  fluctuations  of  ^  is  proportional 

only  to  the  surface  displacement  c(X,Y,t).  The  total  phase  fluctua¬ 
tion  between  rays  reflected  at  surface  level  c(X,Y,t)  and  rays 

reflected  from  the  flat  surface  plane,  i.e.  C  *  0,  is 
0) 

g  ■  -  2  sin  C(X,Y,t).  Figure  2.3  is  a  plot  of  the  two 

£  C  U 

reflected  rays.  If  c(X,Y,t)  is  replaced  by  its  r.m.s.  value  a 

then  gff,  written  as  g,  is  called  the  Rayleigh  parameter,  and  is  given 

by 

4ir a  sin  ij»n 

g- - - -  (2.2.9) 

X 


where  X  =  27Tc/u!q  is  the  acoustic  wavelength.  The  Rayleigh  para¬ 
meter  [75]  is  commonly  used  as  a  measure  of  surface  roughness  ;  g  «  1 
is  a  relatively  smooth  surface  while  g  »  1  corresponds  to  a  rough 
surface. 

The  product  term  — ~  in  equation  (2.1.6)  can  be  expanded 

r0  1 

in  a  similar  way.  From  equations  (2.2.3)  and  (2.2.4),  with  n  being 
set  to  -  *s  ,  we  have 


_ 1_ 

Vl 


1 


)X  +( 


(2.2.10) 
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The  factor  -  needs  to  be  expanded  only  to  first  order  in 

Vi 

X  and  Y  since  it  affects  only  the  magnitude  of  the  integrand  of 

equation  (2.1.6)  whereas  the  term  r^-hr^  determines  the  phase.  Small 

changes  in  r^-H:^  can  produce  large  changes  in  phase,  but  small 

changes  of  — - —  only  cause  a  relatively  small  perturbation  of 

r0rl 

the  integrand.  Also,  because  of  Inequality  in  equation  (2.2.5), 
the  dependence  of  — —  on  £  is  negligible. 

ro  1 

We  find  it  convenient,  furthermore,  to  treat  equation  (2.2.10) 

as  the  first-order  expansion  of  an  exponential  function  so  that 

another  approximation  of  is 

*0  1 


_1 _ 

Vl 


— V  " exp  { ( 
00  10  1 


*00 


*R 


10 


)x  +  ( 


r>* 


(2.2.11) 


’bo 


10 


Substituting  equation  (2.2.1)  into  equations  (2.2.7)  and 


(2.2.11),  rQ  +  rx  and 


can  be  rewritten  in  following  forms: 


0  1 


r^-h^  -  ^Q^IO  +2  ?  Sln  ^0  +  ^[x2  ^1-cos2^q  c°s24>r)-+-Y2  Ci-cos2^  sin2<t>R)] 

(2.2.12) 

-  ~z —  exp  i  - (  ~r - - -)[X  cos  iji  cos  (j>  -Y  cos  i|i sin  <(>_]> 

Vl  r00r10  *  r00  r10  0  R  OR) 

(2.2.13) 


where 


R  - 


2r00r10 

roo+rio 


(2.2.14) 
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2.3  Beam  pattern  function 

The  acoustic  source  is  assumed  to  be  circular-piston 
pressure  transducer  of  diameter  D^.  The  beam  pattern  function 
for  this  type  of  transducer  can  be  written  as  [72,  p.  102] 


Bt(9) 


2J1(^~  Sin  0  ) 


W0^T 

2c 


sin  6 


(2.3.1) 


where  0  is  the  angle  measured  from  the  center  axis  of  the 
piston.  J^( •)  is  a  first  order  Bessel  function,  having  the 
power  series  expansion 


JL(x) 


„  2n 

w  (  —  )  (-1)” 

X  £  *  2  '  * 

2  n-0  (n+1) In! 


(-) 

f  d--^- 


<!> 


312! 


(2.3.2) 


Hence  BT(0)  can  be  rewritten  as 


“0DT  2  W0DT 

(  Bin  0)  <-r~  sin  0) 

Bt(9)  -  1  -  — ^ - +  _2c - 


192 


(2.3.3) 
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If  we  assume  0  is  sufficiently  small  such  that  sin  0  ■  0, 
and  the  beam  function  can  be  approximated  by  the  first  two  terms 
of  equation  (2.3.3),  then 


Bt(6)  -  1  - 


2  2  2 
***0  DT  8 

16c2 


2  2.2 
%  °T  9 

16c2 


(2.3.4) 


Ue  assume  now  that  the  projection  of  the  beam  pattern  on  the 


surface  can  be  approximated  by  a  two-dimensional  Gaussian  function 


Bt(u,v) 


42V 


/  2^.2  .  2.  . 

(u  +v  sin  4»t) 


(2.3.5) 


where  is  the  grazing  angle,  and  where  u  and  v  are  transformed 
coordinates  given  by 


u  ■  -X  cos  -  Ysin  (2.3.6a) 

v  -  -X  sin  +  Ycos  $T  (2.3.6b) 


We  let  the  angle  0  measured  along  the  u  and  v  directions  be  6u 
and  0v  respectively.  Then  for  small  0,  equation  (2.3.5)  can  be 
rewritten  as 


00 


Bt(u,v)  -  e 


lt2(V 


(0  2  +  0  2  sin2  *_) 
u  v  TT 


(2.3.7) 
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By  comparing  equations  (2.3.4)  and  (2.3.7)  and  assuming  a  circular 
beam  pattern,  we  see  that 


2  2 
*0  PT 

16c2 


00 


lt2 (wq) 


or 


W 


4cr 


00 


2c  r 


00 


“oDT 


"Vt 


(2.3.8) 


(2.3.9) 


Substituting  this  into  equation  (2.3.5),  we  have 


Bt(u,v)  »  exp 


_2f  2 

n  fo 


4c 


T  .  2  .  2  .  2.  > 

- 2  (“  +  v  sin  4»t> 

roo 


(2.3.10) 


By  use  of  equation  (2.3.6)  the  beam  pattern  function  can  be  expressed 
in  terms  of  X  and  Y.  Using  Bs  to  indicate  the  source  beam  pattern 
function  on  the  surface,  we  have 


„2f  2  D  2 

B  (X,Y)  -  expj - [X2  (l-co824»Tsin2<|>T)+Y2  (1-cos24»tcos2(|>t) 

4c  rQ0 


+  XY(cos  ^Tsin24>T)  ] 


(2.3.11) 


A  similar  expression  is  obtained  for  the  receiver  beam  pattern  function 


w2f„2  D,  2 


Br(X,Y)  -  exp  |  ~  ^  2  *  2  fx^^“cos^R8in2^R)+Y20-“c°82'|'Rcos2^R) 


‘10 


+  XY(cos  K/Rsin2$R)]  (2.3.12) 


Thus  the  form  of  the  combined  beam  pattern  response  of  both  source 
and  receiver  is 


BsBR  "  eXp  -  —2 


X _ Y_ _ 2XY 

Lxx  Sy  ha  ^ o * 


(2.3.13) 


where 


^  (“0)  2  2  1  2.  2  2.  .  +  2.  2  ,  2  J  (2.3.14a) 

it  fQ  Dt  (1-cos  ^sin  (f>T)  Dr  (1-cos  i|»Rsin  <|>R) 


lyy‘  («n)  =  -5— sr  [  —5 - 5s* - 5 - +  —5 - ^ - 5 - ]  (2.3.14b) 

*  f0  dt  (l-cos^cos  (|>t)  DrZ(1-cosZ*rcos\) 


tJ 2,  ,  ,  r00  +  .  „11Q _ 

°  n2f02  DT2cos2<|»Tsin24>T  DR2cos2i|)Rsin2<|>R 


(2.3.14c) 


For  the  simplified  situation  where  the  source  and  the  receiver 


are  in  a  forward  specular  reflection  geometry,  i.e.  »  $R  ■  » 

and  ■  $ D  +  180°,  equation  (2,3.14)  can  be  simplified  as 

1  K 


» 


(2.3.15a) 


2  l2(“0} 

ha 


2  2 

1* cos  I|»0sin  <|>R 


l2<“0) 


2  2 

l-cos  iJ)qCos  $r 


(2.3.15b) 


ha 


L2(»0) 


2  2 
cos  lpQSin  2«|)r 


(2.3.15c) 


where 


l2<u0> 


(2.3.16) 


Ue  now  substitute  equation  (2.2.12),  (2.2.13)  and  (2.3*13) 
into  equation  (2.1.6)  and  convert  the  integration  over  s  to  one  over 
X,Y  by  using  the  transformations  (see  appendix  B) 


ds 


l+(  IS.  )  +  (  Ji  ) 

A K  3X  1  y  3Y  1 


dXdY 


(2.3.17) 


and 


3_ 

3n 


[i+<  H  )2+  <  If  )2]  1 


ii  1_ 

3X  3X 


ii  i_  i 

3Y  3Y  1 


(2.3.18) 


Then  equation  (2.1.6)  takes  the  form 
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c 


H(ojQ,t)* 


I  /  l3c  3X  9X  3Y  3Yje 
r^r,„  J  J 


7  7  ^fl 

-AX-BY-CXY-DX-EY4-j^F5 


00  10 


where 


2 _ 2 . 


A  "  -J  (1‘cos  VOS  V  +  :  2 


LXX  ^w0) 


B  -  - 


j  (l-co.2Vin2*R)  +  — § 


LyY 


2u>„ 


C  -  - 


j  _0  CO82*0flin*RCO8*R+  -_| 


ScY 


D  “  (  ~~  -  ~)  cosi^ces*- 

roo  rio  0  R 


E  -  -  (  ~~  -  — )cos\|i_sin<|> 
r00  r10  0  R 


F  -  2  sin  1>, 


t  . 


roo  +  r10 


dXdY 


(2.3.19) 


(2.3.20a) 


(2.3.20b) 


(2.3.20c) 


(2.3.20d) 


(2.3.20e) 


(2.3.20f) 


(2.3.21) 


where  L^OjJq) ,  L^Cwq)  and  L^^q)  are  defined  in  equations  (2.3.14). 
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Chapter  III.  Deterministic  Surface 


3.0  Introduction 

The  application  of  the  Kirchhoff  integral  to  the  general 
problem  of  surface  scattered  acoustic  signals  was  described  in  the 
previous  chapter.  An  expression  was  developed  in  equation  (2.3.10) 
for  a  general  water  surface.  Actual  ocean  surfaces  are  highly 
random  and  complex,  but  in  this  chapter  we  side-step  the  issue  of 
random  surfaces  by  considering  a  deterministic  one.  (Random  sur¬ 
faces  are  considered  in  Chapter  IV) .  There  is  substantial  precedent 
for  considering  deterministic  surface  models  as  a  way  of  approaching 
a  solution  for  the  random  surface.  For  instance,  as  mentioned  in 
Chapter  I,  the  early  work  of  Rayleigh  [21]  and  the  more  recent  work 
of  Uretsky  [22],  Murphy  and  Lord  [26],  Heaps  [30]  and  Parker  [76] 
used  a  sinusoidal  surface  to  derive  the  surface  reflection  coefficients. 
Meecham  [77]  considered  the  transmission  of  acoustic  wave  through  a 
sinusoidal  surface  by  the  use  of  the  Kirchhoff  integral  method. 

The  experimental  work  of  Gulin  [36],  and  of  LaCasce  and  Tamarkin  [78] 
dealing  with  reflections  from  pressure-release  sinusoidal  surfaces 
shows  that  much  insight  is  gained  by  analyzing  simple  surface  models. 

If  the  surface  scattering  process  can  be  regarded  as  linear  (as  is 
always  assumed  in  this  work) ,  more  complicated  surface  effects  can  be 
obtained  by  Fourier  analysis  of  the  surface  and  superposition.  This 
was,  for  instance,  done  by  Beckmann  [60]  in  a  study  of  backscatter 
from  a  composite  rough  surface.  He  assumed  each  of  the  component 
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waves  to  be  an  independent  random  process,  and  found  that  the 
small  structure  of  the  surface  (the  component  waves  of  small 
amplitude)  played  an  important  role  in  the  determination  of 

j 

backscattered  acoustic  power. 

In  this  chapter,  we  will  demonstrate  the  effect  of  small 
surface  structures  on  the  frequency  spreading  function  by  considering 
a  wind-blown  surface.  A  commonly  observed  characteristic  of  wind¬ 
blown  surfaces  is  that  the  waves  are  steeper  on  the  downwind  side 
than  on  the  upwind  side  [70].  A  typical  wave  shape  is  shown  in 
Figure  1.2.  A  deterministic  model  that  has  the  feature  of  shallower 
slopes  on  one  side  and  steeper  slopes  on  the  other  can  be  constructed 
with  a  sine  wave  plus  its  second  harmonic;  i.e. 

z(t)  -  h^  cos(wt+$)  +  h,,  os  2u»t  (3.0.1) 

Typical  waveforms  generated  from  equation  (3.0.1)  for  different 
values  of  h^,  h^  and  $  are  shown  in  Figure  3.1.  It  is  seen  that  by 
proper  choices  of  amplitudes  and  phase  angles  (e.g.  =  5h2»  4>“45°) , 

this  two-harmonic  surface  wave  will  have  a  slope  distribution  similar 
to  that  of  a  wind  generated  wave. 

The  frequency  spreading  function  of  the  acoustic  signal 
scattered  from  this  deterministic  surface  model  is  discussed  for  two 
scatter  geometries,  crosswind  and  down(up)  wind.  Crosswind  means 
that  the  wind  direction  is  at  right  angles  to  a  vertical  plane 
containing  source  and  receiver,  and  down  (up)  wind  indicates  that 
the  wind  direction  is  parallel  to  this  plane.  Thus  in  Figure  3.2, 

0°,90°(270°) ,  and  180°  correspond  respectively  to  upwind,  crosswind 
and  downwind. 
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3.1  The  surface  model 


( 

In  accordance  with  our  discussion  in  section  3.0,  the  surface 
displacement  at  an  arbitrary  point  of  a  wind  generated  surface  is 
described  by  a  time  function  as  in  equation  (3.0.1).  If  we  choose 
our  coordinate  system  such  that  the  wind  blows  in  the  positive  Y 
direction  and  if  the  surface  wave  is  one-dimensional,  the  surface 
wave  motion  can  be  expressed  in  the  following  form 


C(X,Y,t)  =  h-jCOsfpY-ft-H)  +  h2cos(2pY-20t) 


C 

P 


_ft 

P 


(3.1.1) 

(3.1.2) 


where  p  is  the  surface  wave  number,  ft  is  the  wave  frequency,  and 
Cp  is  the  phase  velocity. 

Observe  that  the  fixed  phase  angle  $  in  equation  (3.1.1)  and 

the  constant  phase  velocity  for  both  harmonics  is  not  consistent 

with  the  model  usually  assumed  for  gravity  waves  (e.g.  Parkins  [A] 

«L2 

or  Scharf  and  Swartz  [5]).  For  gravity  waves,  p^  *  — — —  ,  and  the 
phase  angle  is  usually  regarded  as  being  uniformly  distributed 
between  0  and  2ir.  The  dispersion  relation  for  gravity  waves  does  not 
permit  persistent  slope  asymmetries  of  the  kind  required  in  our  dis¬ 
cussion  here.  Equation  (3.1.1)  therefore  should  be  regarded  as  a 
rather  idealized  model  for  a  wind-blown  surface. 

In  order  to  study  the  surface  scattering  behavior  from  this 


model,  it  is  necessary  to  replace  s(X,y,t)  of  equation  (2.3.19) 
with  the  expression  in  equation  (3.1.1)  and  then  evaluate  the 
integral.  We  first  consider  the  downwind  case. 


1 


3 . 2  Downwind 


Figure  3.3  Illustrates  a  downwind  *  180°)  scatter  geometry. 
The  acoustic  source  is  on  the  -Y  side  and  the  receiver  is  at  the 
+Y  side.  As  we  showed  in  Chapter  II,  the  distances  from  an  arbitrary 
surface  point  to  the  source  and  receiver  are  respectively  r^  and  r^. 
These  have  the  forms: 


r0  -  ft  +  <Y«00cos*0)2  +  tt«00.in*0)2  (3.2.1) 

-  ft  +  <Y-r10cos+0)2  +  (£+r10sini|.0>2  (3.2.2) 


By  the  use  of  the  Fresnel  approximation  as  described  in 
Chapter  II,  the  total  transmission  distance  r^  +  r^  can  be  approximated 
as  equation  (3.2.3).  From  equation  (2.2.7),  we  have 

rQ  +  rx  -  r0Q  +  r1Q  +  2^  sim|»0  +  ^  [X2+Y2sin2<|>0]  (3.2.3) 


From  equation  (2.2.11),  the  approximation  of  the  attenuation  factor 
(r^r^)-^  can  be  put  into  form: 


— i  -  — -  exp  -I  -  (  — - - )Y  cos  i|»0j 

r0rl  roorio  r00  r10 


(3.2.4) 
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Both  equation  (3.2.3)  and  (3.2,4)  are  good  approximations  if  the 
following  conditions  hold: 


C(X, Y,t)  «  X,  Y  (3.2.5a) 

X»Y  «  r^,  r1(J  (3.2.5b) 

In  the  mathematical  preliminaries,  we  showed  that  the 
projection  of  the  combined  source-receiver  beam  response  on  the 
X-Y  plane  can  be  approximated  by  a  Gaussian  function.  For  a 
downwind  geometry,  this  combined  beam  response  can  be  expressed 
in  a  simple  form  as  equation  (3.2.6)  by  replacing  <j>R  with  180° 
in  equation  (2.3.15) 


B  B  ■  exp(  -  (X2+Y2  sin2*  )j  (3.2.6) 

1  LZ(u>0)  0  » 

Thus  L(Uq)  is  the  acoustic  beamwidth  measured  in  the  X  direction. 

In  order  to  compute  the  power  spectrum,  we  first  have  to 
evaluate  the  transfer  function  H(u>Q,t)  in  equation  (2.3.19).  After 
substituting  equation  (3.2.3,  4,  6)  into  equation  (2.3.19),  H(u^,t) 
can  be  rewritten  as 


H(u0,t) 


jVs 


4,rroorio 


~ — (x2+y 

L  <V 


a)  X2+Y2sin2* 

-  (  — - ^Ycos*  +  j  [ - 

r00  r10  R 


-  +  2C  sin*0] |  dXdY . 


sin 


In  this  section  we  assume  that  the  surface  height  C(X,Y,t)  is 

small  compared  with  the  incident  acoustic  wavelength  so  that 
“o 

2  —  s(X,Y,t)sin  iJ/q  <<  1.  Thus  the  phase  term  can  be  approximated 
by  its  linear  expansion,  i.e„ 


j2  C(X,Y,t)sin<|/  u 

e  c  0  -  1  +  j2  ^  C  (X,Y,t)sin 


(3.2.8) 


Observe  that  the  surface  wave  motion  is  assumed  to  be  one- 
dimensional  and  is  propagating  in  the  Y  direction.  Therefore,  the 
slope  is  zero.  We  replace  £(X,Y,t)  by  the  expression  in  equation 
(3.1.1)  and  substitute  equation  (3.2.8)  into  H(o>q,  t) .  The  integration 
of  equation  (3.2.7)  is  straightforward  but  lengthy.  Details  of  the 
computation  are  contained  in  Appendix  C.  The  result  of  the  integration 
can  be  put  into  the  form: 


H(wQ,t) 


(3.2.9) 


Expressions  for  A,  B,  Hq  and  are  also  contained  in  Appendix  C. 
The  expression  of  the  received  signal  r(t)  is  the  product 

-jv 

of  H(Ug,t)  and  e 


-JV 

r(t)  -  H(u>0,t)e 


-jwnt  J5,  - j  (wn+mft)  t  -j(u>n-m£l)t 

-  do  •  +  £iiv •  +  ve  ’  a-2-10> 
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Equation  (3.2.10)  indicates  that  the  interaction  of  an  acoustic 
sinusoidal  wave  of  frequency  with  a  moving  surface  corrugation 
given  by  the  expression  in  equation  (3.1.1)  results  in  the  generation 
of  frequencies  u)q  +  mfi  with  m  ranging  from  zero  to  at  least  6. 

The  number  6  arises  from  the  first  order  expansion  of  the  phase  term 

“o 

exp  [-j2  —  sin  1^5]  .  For  surfaces  of  greater  roughness,  a  higher 
order  expansion  of  the  phase  term  would  become  necessary,  and  would 
result  in  more  frequency  components. 

From  equation  (3.2.10)  we  see  that  the  coefficient  Dq  Is 
associated  with  the  transmitted  frequency  Wq  and  represents  the 
amplitude  of  the  coherent  (or  specular)  surface  reflection.  If  we 
neglect  the  beam  pattern  effect  in  A  and  B,  l.e.  if  the  beam  width 
L(uq)  in  equation  (c“2)  is  assumed  to  be  infinity,  then  Dq  can  be 
written  as 


-ju>n(t-T  ) 


r00+r10 


2.  2 

h 

2 


3  u)r 


8  l  P  n  r\  n  A 

[1 - 2p  n,  ]-j - -  p  hx  h2sin  ipQcos  2<J>! 


’  exp|\ 


2  c 

(rio-roo>2  cot2  * 


2  c  r00r10<t10+r00) 


(3.2.11) 


The  above  equation  is  obtained  from  equation  (C-ll) .  The 
leading  term  on  the  right-hand  side  of  equation  (3.2.11)  represents 
the  phase  change  and  attenuation  due  to  the  transmission  distance 
roo  +  r-in  1°  the  specular  reflection  path.  The  term  in  the  square 


l10 


brackets,  1  - 


A  2 

P  hl 


.  2.  2 


“  2p  h2  ,  is  an  additional  attenuation  of  the 


3-11 


smooth  surface  response  due  to  the  surface  corrugation.  The  factors 


ph^  and  are  proportional  to  the  r.o.s.  slope  of  the  first  and 

second  harmonic  of  the  surface  wave.  It  is  reasonable  that  the 

larger  the  corrugation  is,  i.e.  larger  ph^  and  ph^,  the  more 

the  coherent  component  of  the  acoustic  signal  is  attenuated.  The 
3  w0  2  2 

Imaginary  term  j  j  —  P  ^  h2Sin  t^cos  2$  is  a  phase  lag  which 
is  a  function  of  the  phase  angle  $.  It  is  zero  for  4  “  45°  and 
is  maximum  for  $  ■  0°.  From  the  surface  profile  shown  in  Figure 
3.1,  we  observe  that  $  ■  45°  results  in  a  "zero  mean"  sawtooth 
wave  from  which  one  would  expect  no  net  phase  shift.  However, for 
$  ■  0°  (or  180°) ,  the  surface  wave  is  characterized  by  a  sequence 
of  narrow  pulses,  and  the  average  height  of  the  reflecting  facets 


can  therefore  be  above  or  below  the  smooth-surface  reference.  The 

r  <ri<fr00)2  cot2*«  1 


last  phase  term  exp  £ j 


2  c  roorio(rio+roo) 


represents  the 


phase  change  due  to  scatter  geometry.  If  the  scatter  geometry  is 


symmetric,  i.e.  r^  -  r^,  then  the  argument  of  this  phase  term 
becomes  zero.  This  implies  that  the  coherent  reflection  is  associated 


with  the  specular  ray-path.  When  the  scatter  geometry  is  asymmetric, 
i.e.  r^Q  \  r^Q,  then  there  is  a  net  phase  shift  in  addition  to  the 
nominal  phase  change  UqTs  .  Hence  under  this  condition  coherent 
reflection  is  no  longer  associated  with  the  specular  ray.  This 
effect  can  also  be  explained  by  considering  the  Fresnel  zones  [75]. 

A  symmetric  downwind  scatter  geometry  is  shown  in  Figure  3.4.  The 
ellipses  on  the  surface  plane  represent  the  Fresnel  zones.  Figure  3.5 

is  the  similar  plot  as  Figure  3.4  except  the  scattering  geometry  is 


ia  asymmetric  with  .  In  Figure  3.4  the  specular  point 

la  at  the  center  of  the  first  Fresnel  zone  (the  inner-most  ellipse) 
and  therefore  the  specular  path  is  associated  with  the  coherent 
reflection.  In  Figure  3.5,  the  center  of  the  first  Fresnel  zone 
is  shifted  away  from  the  specular  point  and  toward  the  receiver 
(toward  the  source  if  r^^  >  r  ) .  Thus  there  is  a  net  phase 
difference  between  the  coherent  reflection  and  specular  reflection. 


The  coefficients  of  the  side-frequencies  (the  frequency 
component  other  than  the  transmitted  frequency)  above  and  below 
are  respectively  and  Dm_,  m  «  1,6.  Examination  of  these 
coefficients  [cf.  equation  (C-12)]  reveals  that  these  coefficients 
with  m  ranging  from  1  to  4,  are  in  general  not  complex  conjugates  of 
each  other.  Therefore,  the  magnitudes  of  D  .  and  D  are  not  equal 

nfr  in— 

for  m  ■  1  to  4.  We  are  mainly  interested  in  the  power  of  the  first- 

^  2  2 

order  side-frequencies  ,  i.e.  |D^+|  and  | D^_[  .  These  can  be 
written  in  the  forms: 


lDi+(Vt)|2  “  [l2  +  **+X-T  sin  24 J  *  |  exp 
|Dl_(ci)Q,t)|2  -  X2  +  Y2-  J  Y  sin  24J'  |exp[ 


(E-jp)2 

2B 

(E+jp)2 

2B 


]l 

]l 


(3.2.12) 


(3.2.13) 


where 


'oo^io 


2  h.  «  Z 

(1-p  (  -^-  +  2h* - § - 5—)] 

4  4mq  ala 


(3.2.14a) 


Because  of  the  first-order  expansion  of  equation  (3.2.8),  the  higher- 
order  side-frequency  coefficients  are  probably  not  correct  in  any 
case. 
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igure  3.5  Fresnel  zones-asymmetric  source-receiver  geometry  (r 


y  5  p2hlh2 

"  4  '  roo+Tio 


(3.2.14b) 


and  where  from  the  definition  In  appendix  C,  l.e.  equation  (C- 2) 


u 


B  -  [ 


l2<»0) 


■151  sl“2*o 


(3.2.15a) 


E "  ( ^r~  ■  )  c°8  *0 

00  10  u 


(3.2.15b) 


2  2 

Notice  that  |D^+(a>Q,t)  |  and  |Dj_(a>Q,t)|  are  in  general  not 
equal.  They  differ  because  of  the  sign  differences  in  the  bracketed 
term  and  in  the  exponential  term.  These  two  sign  differences 
represent  two  different  mechanisms.  The  sign  inside  the  brackets, 
associated  with  the  term  X  T  sin  2$ ,  is  the  result  of  asymmetric 
water  surface  slopes  (the  assumed  wind-driven  surface  model) .  The 
second  sign  difference  in  the  exponential  term  results  from  the 
asymmetric  scattering  geometry  (i.e.  r^  +  r^)  . 

We  consider  first  the  case  rnn  ■  r,n.  This  results  in 

uu  xu 

E  *  0  [cf.  equation  (3.2.15b)]  and  leaves  only  the  sign  difference 
in  the  bracketed  term.  We  observe  that  the  asymmetry  (the  amplitude 
difference  between  |Dj+(u)Q,t)|  and  |Dj_(<i>Q,t)|  exists  when  $  ^  0 
(or  ♦  ■  ip-  ,  n  ■  integer),  and  is  most  pronounced  when  ♦  -  45°. 

In  Figure  3.1,  $  -  45°  represents  the  most  asymmetric  saw  tooth 
wave.  There  is  no  asymmetry  when  the  amplitude  of  the  second  harmonic 
is  zero  (i.e.  h^  ■  0,  T  «  0) .  The  question  as  to  which  one  of  the  two 
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side-frequencies  has  the  larger  power  magnitude  depends  on  the 
sign  of  X  +  .  Generally  speaking,  the  water  surface  under 
consideration  has  no  surface  slope  with  magnitude  larger  than  1. 
This  condition  is  required  by  the  laws  of  hydrodynamics  for  non¬ 
breaking  surface  waves.  Hence  7  has  a  positve  magnitude  at  all 
times,  (see  appendix  D.) 

We  can  conclude  therefore  that  under  the  condition  r^  =  r^Q> 

2 

the  upper  side-frequency  power  magnitude  |D.j+(<j>Q,t)  |  is  always 

2 

larger  than  the  lower  side-frequency  power  magnitude  |D^_(u>Q,t)| 
for  0  <  $  <  90°.  i.e.  for  wind  driven  waves. 


Observe  that  the  effect  considered  here  is  a  function  of  the 
surface  slope  asymmetry  only  and  does  not  depend  on  source-receiver 
geometry.  Exchanging  the  positions  of  source  and  receiver  replaces 
$0  by  (see  Figure  3.2),  and  since  tl>Q  enters  the  expression 

only  through  sin  ,  replacing  by  ir-^Q  has  no  effect  on  the 
magnitudes  of  |Dj+(u»Q,t)|  and  |Dj_(u>Q,t)  |  .  This  result  can  also 
be  justified  heurlstically  by  invoking  simple  reciprocity. 


The  power  ratio  between  the  upper  side-frequency  and  the 
lower  side- frequency  can  be  obtained  from  equation  (3,2.12)  and 


(3.2.13) 


Dj+O^.t) 

Dj-O^.t) 


2 


X2  +  Y2  +  2X7  sin  2$ 
X2  +  T2  -  T  sin  2$ 


(3.2.16) 


The  values  of  lie  between  0°  and  90°  in  order  to  have  the  wind 
driven  waves  moving  in  the  direction  of  the  steeper  wavefront. 1 70 J 
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A  plot  of  this  power  ratio  as  function  of  $  according  to 
equation  (3.2.16)  is  shown  in  Figure  3.6. 

We  now  focus  our  attention  on  the  sign  difference  in  the 
exponential  term.  To  this  end*  we  set  $  *  0,  l.e.  to  have  a 
surface  wave  with  symmetric  slope  distribution.  By  replacing 
B  and  E  with  the  expressions  in  equation  (3.2.15)*  the  exponential 
term  becomes  „ 


1  r COS  *0  2]  “o  2p  COS  ^0 


r  (E+jp)2^ 
exp  I 

2  2  "  P 
j-  L2(io0)L  R'2 

1 

R' 

1 

L  2B  J 

■  exp 

2  sin2  T  -A- — 

-+  “°2  1 

J 

where 


(3.2.17a) 


R'  roo  rio 


(3.2.17b) 


The  power  ratio  of  the  upper  and  the  lower  side-frequencies 
can  be  written  as 


D1+(<vt) 

D1.(a>0,t) 


p  cos 


^  2  r  _i _  Mo21 

sin  i|>0  [  4  +  2  2  J 

L  L  (oO  c  R 


(3.2.18) 


Figure  3.7  shows  a  log-log  plot  of  |Dj+(w0,t)/Dj-(u)0,t)  | 
as  a  function  of  rig/rQQ»  with  grazing  angle  as  a  parameter. 

Observe  that  the  power  ratio  is  unity  when  r^  -  r^.  It  is  greater 
than  unity  when  r^  <  r^  .  Also  notice  that  the  power  ratio 
saturates  at  very  large  or  very  small  values  of  ri(/r00*  The  limiting 
value  can  be  obtained  from  equation  (3.2.18)  by  letting  ri(/r00  8° 
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to  infinity*, 

2  A 

lint  VV;  .  r  2t>)QcpL  (at^cos  »Q  (3<2>19) 

^  “  Di-(  0*°  **  8±n2  '<'o{4c2roo2+uo2L4<“o)1 

where  R  and  R'  are  approximated  by  and  r^,  respectively. 

As  far  as  the  effect  of  grazing  angle  ijig  is  concerned,  we  see 

that  for  4*q  “  90°  (vertical  incidence  and  reflection)  the  power 

ratio  is  unity,  but  for  very  shallow  grazing  angles  this  ratio  can 

become  very  large.  Figure  3.8  has  a  semi- log  plot  of 

2 

log | Dj+ (o>q , t) /D^_ (u)q , t)  |  versus  grazing  angle  In  fact,  in 

the  limit  of  *  0°  or  =  180°,  the  power  ratio  goes  to  infinity. 
However,  for  extremely  low  grazing  angles  as  well  as  for  extreme 
values  of  r^/r^,  the  assumption  of  no  shadowing  is  violated,  and 
therefore  the  results  should  not  be  carried  out  to  these  limits. 

The  power  ratio  in  equation  (3.2.18)  is  sensitive  to  the 
wind  direction,  and  will  change  to  its  reciprocal  value  when  the 
direction  of  wind  reverses.  For  example,  if  we  take  the  scatter 
geometry  in  Figure  3,3  with  r^  <  ,  then  according  to  equation 

(3.2.18)  the  power  of  the  upper  side-frequency  is  larger  than  that 
of  the  lower  one.  If  we  reverse  the  wind  direction  so  that  the  wind 
blows  along  the  -y  axis,  we  are  effectively  replacing  the  wave  number 
p  with  -p  in  equation  (3.2.18).  Thus  the  power  ratio  is  inverted 
and  the  lower  side-frequency  will  become  the  one  with  a  larger  power 

*  The  Fresnel  approximation  breaks  down  before  reaching  this  limit. 
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Figure  3. 


First-order  side  frequency  power  ratio  versus 
grazing  angle  <Pq;  the  dashed  curve  represents 
the  power  ratio  due  to  asymmetric  surface  slopes 
for  $  ■  45°. 

(cf.  Figure  3.6) 


magnitude.  This  same  effect  can  also  be  observed  from  Figure  3.3 
and  equation  (3.2.18)  by  letting  exceed  90° ,  l.e.  the  source 
moves  to  the  +y  side  and  receiver  to  the  -y  side,  and  the  wind 
blows  along  +y  direction. 

If  the  beam  width  L(<i>q)  of  the  combined  source-receiver 
beam  response  becomes  very  large,  i.e.  if  both  source  and  receiver 
lose  their  directivities,  and  for  Tqq  <<  r^,  equation  (3.2.19) 
can  be  further  simplified  to  the  form: 


lim 

rl(/r00 

L 


Dj+(w0,t) 

Dl_  (<»»0.  t) 


2  cos  i|i 


exp 


uQ/c  1-cob  i|>q 


H 


(3.2.20) 


We  see  that  the  argument  of  the  exponential  function  consists 
of  two  parts,  one  is  the  ratio  of  the  surface  wave  number  p  to  the 
acoustic  wave  number  Uq/c,  the  other  is  a  function  of  the  grazing  angle 
a  geometric  factor.  According  to  equation  (3.2.20),  the  power 
ratio  is  large  when  pc/ui^  is  large;  and  the  ratio  approaches  unity 
when  pc/(0q  closes  to  zero.  A  heuristic  explanation  of  this  is  that 
for  a  given  acoustic  wave  number,  if  pc/m^  decreases,  the  surface  wave¬ 
length  must  increase.  In  the  limit  pc/oJq  -*  0,  the  surface  wavelength 
is  infinity  and  the  surface  wave  becomes  a  plane  surface.  All 
reflections  become  specular,  and  there  is  no  asymmetry. 

The  effects  of  the  geometry  asymmetry  and  the  surface  slope 


asymmetry  on  the  power  ratio  are  not  equal  in  general.  The  logarithmic 
value  of  the  power  ratio  due  to  the  surface  slope  effect  is  plotted 


> 


dashed  in  Figures  3.7  and  3.8.  We  see  that  the  effect  of  surface 
slope  is  smaller  than  that  of  the  geometry  asymmetry  and  is 
significant  only  when  r^  and  are  approximately  equal  or  when 
^  the  grazing  angle  is  close  to  90°. 

I 

» 

I 

ft 
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3.3  Crosswlnd 


The  discussion  for  crosswlnd  Is  very  similar  to  that  for 
downwind.  The  scattering  geometry  is  shown  in  Figure  3.9  with 
the  source  and  receiver  located  in  the  x-z  plane.  The  mathematical 
operations  are  very  similar  to  those  already  performed  for  downwind. 
Details  of  the  computation  are  given  in  appendix  E.  We  find  that 
the  received  signal  can  again  be  expressed  as  in  equation  (3.2.10); 
however,  the  expressions  for  the  D's  are  generally  different. 

The  coherent  response  Dq  has  the  same  form  as  in  the  previous 
section  and  we  need  not  consider  it  again. 


Expressions  for  power  magnitudes  of  the  first  order 
side-frequencies  are 


|Dj+(ti)Q,t)  |  2  -  {  X2  +  Y2  +  Xf  sin  2$  }  •  exp  { 


„2 

E  -P 


2B 


E2-p2 

iDj-CoiQ.t)  |2  -  {  X2  +  Y2  ~  XY  sin  2*}  •  exp  {  — ^ - 


1  (3.3.1) 


}  (3.3.2) 


where  X,  Y,  B  and  E  are  defined  in  equations  (3.2.14)  and  (3.2.15). 

The  sign  difference  in  these  two  expressions  is  related  to 
the  surface  wave  slope  asymmetry.  The  sign  difference  arising  from 
asymmetry  in  the  scatter  geometry  is  no  longer  present.  Therefore, 
for  a  symmetric  surface,  e.g.,  gravity  waves,  we  would  expect 
symmetric  side-frequencies  in  the  crosswind  geometry  for  any  r^Q 
and  r^Q,  (r^  and  r^  are  assumed  to  satisfy  the  far-fleld  assumptions). 


2353338 
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The  detailed  consideration  of  the  sign  difference  in 
equations  (3.3.1)  and  (3.3.2)  is  the  same  as  discussed  in  the 
previous  section  since  only  the  slope  effect  is  presented  here, 
the  power  magnitude  of  the  upper  side-frequency  is  always  larger 
than  that  of  the  lower  side-frequency.  We  see  that  there  is  no 
effect  of  geometry  asymmetry  in  the  crosswind  condition;  therefore 
this  is  a  convenient  experimental  conflguraton  for  the  observation 
of  surface  slope  asymmetries  [3]. 
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3.4  Summary 

In  this  Chapter  we  have  considered  the  frequency  spreading 
function  for  both  the  downwind  and  the  crosswind  conditions  using 
a  deterministic  model.  We  have  found  that  the  power  magnitudes  of 
the  upper  and  the  lower  side-frequencies  produced  by  this  model  are 
generally  unequal.  There  are  two  different  mechanisms  that  cause 
inequality.  One  is  the  scattering  geometric  asymmetry,  l.e.  the 
source  and  receiver  are  at  different  depths.  If  the  source  is 
located  at  a.  smaller  depth  than  the  receiver  and  if  the  wind  blows 
from  the  source  to  the  receiver  (l.e.  surface  wave  propagates  from 
source  to  receiver) ,  then  the  power  of  the  upper  side-frequency  is 
larger.  This  effect  is  reversed  when  wind  direction  (the  surface 
wave  direction)  is  reversed  or  if  the  wind  direction  is  unchanged 
but  the  source  is  at  the  larger  depth.  This  mechanism  does  not 
exist  in  the  crosswind  condition.  The  other  mechanism  is  the 
surface  asymmetry  caused  by  asymmetric  surface  slopes  in  wind- 
driven  waves.  This,  taken  by  Itself,  always  results  in  a  larger 
power  magnitude  in  the  upper  side-freqenncy .  The  power  ratio 
(asymmetry)  caused  by  asummetric  surface  slopes  is  generally  smaller 
than  that  caused  by  the  asymmetric  scatter  geometry  as  observed  in 
Figures  3.7  and  3.8.  Therefore,  the  overall  power  difference  tends 
to  be  mainly  a  function  of  the  scattering  geometry. 


» 
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Chapter  IV  Random  Surface 


4.0  Introduction 


In  this  Chapter  we  consider  a  random  water  surface  with  a 
Gaussian  surface-height  distribution.  As  already  discussed  in 
Chapter  I,  the  Gaussian  assumption  is  frequently  used  in  the 
analysis  of  surface  scatter  even  if  it  is  not  rigorously  correct 
in  all  cases.  Parkins  [4],  Scharf  and  Swarts  [5],  McDonald  and 
Schultheiss  [10]  used  this  assumption  to  study  the  frequency 
spreading  in  the  forward  acoustic  scattering  from  random  rough 
surfaces.  Schwarze  [80]  developed  the  doppler  spread  of  back- 
scattered  acoustic  wave  from  e  Gaussian  sea  surface.  Clay  and 
Medwin  [7]  computed  the  covariance  of  signals  observed  at  a  pair 
of  receivers  with  the  aid  of  Gaussian  surface.  Beckmann  [49] 
has  considered  the  surface  scattering  by  non-Gausslan  surfaces. 

He  found  that  the  results  in  most  scattering  analyses  are  not  very 
sensitive  to  the  precise  choice  of  distribution. + 

The  analysis  considers  a  general  scattering  geometry  with 
arbitrary  source  and  receiver  locations.  Thus  the  study  includes 
various  oblique  and  asymmetric  scattering  situations.  In  the 
previous  Chapter,  we  saw  that  surface  slopes  played  an  important  role 
in  determining  the  shape  of  the  frequency  spreading  function.  Hence 
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slope  staCisti.cs  are  explicitly  Included  in  this  analysis.  By 
assuming  independence  of  surface  slope  and  surface  height 
statistics,  asymmetry  effects  in  the  received  signal  spectrum 
that  are  similar  to  those  obtained  for  deterministic  surfaces  can 
be  shown  to  exist  also  for  random  surfaces. 


4.1  General  scattering  geometry 


If  the  water  surface  Is  random,  the  surface  scatter  process 
can  be  regarded  as  a  randomly  time-varying  linear  system  having  a 
transfe*  ^unction  H(u)Q,t).  This  represents  the  instantaneous 
amplitude  and  phase  at  the  receiver  due  to  a  sinusoidal  signal  of 
frequency  from  the  source.  The  received  signal  has  the  form 

r(t)  =  H(u>0,t)e  u  (4.1.1) 

Figure  4.1  shows  the  scattering  geometry.  The  origin  of  the 
coordinate  system  is  chosen  to  be  the  aiming  point  of  the  acoustic 
source  on  the  flat  surface.  In  other  words,  the  source  has  a  beam 
pattern  whose  center  line  intercepts  the  flat  surface  at  the  origin. 
The  x,y  and  z  axes  are  as  previously  defined.  The  location  of  the 
acoustic  source  is  specified  through  three  polar  coordinate  para¬ 
meters,  Tqq,  and  .  Tqq  is  the  distance  from  source  to  system 
origin,  <(>T  is  the  azimuthal  angle  of  r^  projection  on  the  x-y 
plane,  measured  positive  clockwise  from  the  -y  axis,  is  the 
grazing  angle.  The  receiver  location  can  be  described  through  a 
similar  set  of  parameters  r^^,  <f>R  and  i|>R. 

The  transfer  function  for  this  scattering  geometry  can  be 
derived  from  the  Fresnel-corrected  Kirchhoff  integral  and  has  the 
form  (see  Appendix  F) 
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General  scatter  geometry  with  wind  direction 


H(u>0,t) 


1  (r00+r10^  f  f  -Ax2 -By 2 -C xy-Dx-Ey+j-^F^ 

~ - - - —  /  /  {  t~  -?  I— 3e  dxdy 

4,  r  r  -«  «  ^  x3x  y3y 

47r  r00r10 

(4.1.2) 


where  A,  B,  C,  D,  E  and  F  are  specified  in  Appendix  F. 

Our  interest  is  in  the  power  spectrum  of  the  received  signal 
due  to  a  sinusoidal  input.  The  power  spectrum  can  be  obtained  by 
Fourier  transforming  the  autocorrelation  function  of  the  received 
signal  r(t).  By  the  assumption  of  wide-sense  stationarity  in  the 
acoustic  scatter  channel,  we  can  write  the  autocorrelation  function 
in  the  form  [82] 


r(t)  r*(t  +  x) 


$(td0,t)e 


(4.1.3) 


where  $((j>q,t)  «  H(uiq, t)H*(a)Q,t  +  x) 


(4.1.4a) 


is  the  autocorrelation  function  of  the  scatter  channel.  Substituting 
the  expression  of  H(iOg,t)  in  equation  (4.1.2)  into  equation  (4.1.4a), 
we  get  a  four-fold  integral  as  follows 


oo  *  oo  -  oo  *  ao 


■  TTr-’V’T  l ill 

loir  Tqq  r^Q  -oo  -a>  -oo  -oo  1  2  12 


(4.1.5) 


O 


where 


J<5V  *?v  } 

X1  2  yl  y2 


{j  —  F(l-C  2-S  2)+  (2Ax1+Cy1+D)+C  (2By.+Cx1+E)  } 

J  c  yx  x1  1  '1  yx  1  1 


■{-i^F(l-C  *2-  K  *2)+?  (2A*x.+C*y  +D*)+C  (2B*y_+C*x  +E*)  } 

c  x2  y2  x2  2  2  y2  1  l 


•  exp{ -Ax12-By 12-Cx1y1-Dx1-Ey 1-A*x24-B*y2i-C*x2y2-D*x2-E*y2  } 

(4.1.5a) 


2  2 


q<c1*c2)  =  exp{j  ~r  F^1<xi»yi»t>- 1  F  c2(x2,y2*t  +  T)  1  (4-1-5b) 

We  use  the  subscripts  1  and  2  to  indicate  variables  at  time 

t  and  t  +  t  respectively.  Notice  that  J(t  )  summarizes 

1  2  yl  y2 

the  channel  slope  information  and  Q^,^)  is  a  function  of  surface 

waveheight.  The  averaging  operation  in  the  integral  is  over  the 

product  of  J (C  ,C  ,C  )  andQC^.C  ).  The  random  variables 

X1  x2  yl  y2  L 

to  be  averaged  are  the  surface  slopes,  t  ,S  ,C  •  and  surface 

X1  x2  yl  y2 

displacements  and  C2« 

The  averaging  operation  is  obviously  fairly  complicated.  In 
order  to  simplify  this  operation  we  make  two  assumptions.  The  first 
/{  these  is  to  assume  a  one-dimensional  random  surface  with  randomness 
n.v  In  the  wind  direction.  In  this  way  the  cross-wind  slope  terms 
«nu  - _  have  negligible  effect  on  the  final  result  and  can  be 

*  *  y 

. « r  . .  This  approximation  can  be  justified  by  the  final  result 
,v  the  effect  of  all  of  the  slope  terms  is,  in  fact, 
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quite  small.  We  do  not  wish  to  eliminate  the  downwind  slope 
terms  since  they  contribute  an  asymmetry  to  the  frequency  spreading 
function,  whose  study  is  one  of  the  objectives  of  this  analysis. 

Even  though  this  is  a  small  effect  it  is  important.  The  crosswind 
slopes  do  not  contribute  such  an  asymmetry  and  setting  them  to  zero 
therefore  does  not  alter  the  final  result  in  any  essential  way. 

The  second  assumption  concerns  the  averaging  operation  itself. 

It  can  be  performed  exactly  only  if  the  surface  slopes  and  heights 
are  jointly  Gaussian.  However,  the  Gaussian  assumption  is  inconsistent 
with  the  directional  asymmetry  of  windblown  surfaces  that  we  believe 
to  be  partly  responsible  for  the  phenomenon  being  studied. here.  This 
follows  from  the  fact  that  a  strictly  Gaussian  surface  can  be 
represented  by  a  two-dimensional  "Rice"  model  [83];  i.e.  an  expansion 
in  terms  of  sinusoids  with  uniformly  distributed  random  phase  angles. 
Such  a  surface  would  be  statistically  isotropic.  Hence  a  more 
fruitful  simplifying  assumption  for  the  evaluation  of  equation  (4.1.5) 
is  that  the  surface  heights  are  approximately  Gaussian  and  that  they 
are  uncorrelated  from  the  slopes.  A  possible  basis  for  the  second 
part  of  this  statement  is  that  the  slopes  depend  largely  on  small 
ripples  of  the  surface,  when  the  surface  is  rough,  and  that  ripples 
of  any  slope  can  occur  at  any  surface  level.  This  may  not  be  quite 
true  for  the  small-surface-roughness  situation.  However,  if  both 
surface  height  and  slopes  have  small  magnitudes,  cross-product  terms 
can  be  regarded  as  second-order  magnitudes.  Hence  the  assumption  of 
uncorrelated  surface  heights  and  slopes  results  in  an  approximation 
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that  is  correct  to  first  order. 
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With  these  two  assumptions,  the  integral  of  equation 

(4.1.5)  can  be  rewritten  as 


*(o>0,t)  =  - ^ - 2 - 2 


00  OO  00  00 

//// 


1611  r00  rlO  -  -  -  J(S1‘S2)‘Q(H,S)dXldX2dyldy2 


(4.1.6) 


where 


J^y  > 

yl  y2 


wo  ?  “n  ? 

(j  ~  F(l-C  )+C  (2By  +Cx  +E)}*{-j  ~r  F(l-C  *  )+Z  *(2B*y.+c*x  +E*)j 

yl  yl  c  y2  y2  £‘ 


•exp{ -Ax12-By12-Cx1y1-Dx1-Ey1-A*x22-B*y22-C*x2y2-D*x2-E*y2 } 

(4.1.7) 

The  integrand  is  the  product  of  two  averaged  terms.  One  is 

J(C  , ?  ),  which  contains  the  surface  moments  up  to  the  fourth 

yl  y2  _ 

order.  The  other  is  Q(^,^2),  which  is  the  characteristic  function 
of  the  Gaussian  distributed  variables  4^  and  42»  It  can  be  written 
in  the  form 

_  2 

Q(41,C2)  -  exp{-  -y-  F2a2[l-y(5,n,T)]}  (4.1.8) 

c 

where  a  is  the  r.m.s.  surface  waveheight  and  Y(£>n,T)  is  the 
normalized  surface  height  correlation  function,  given  by 

^(C.n.T)  -  51(x1,y1,t)C2(x2»y2't+T^  (4.1.9) 
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where 


(4.1.10a) 


5  -  x2-x2 


n  -  y2-yx 


(4.1.10b) 


The  surface  correlation  function  Is  a  critical  factor  In  the 
final  result.  The  expression  chosen  for  It  should  represent  the 
actual  surface  with  some  degree  of  realism,  and  at  the  same  time 
It  should  be  simple  enough  to  permit  evaluation  of  equation  (4.1.6). 
The  correlation  function  of  the  ocean  surface  has  been  measured  by 
a  number  of  researchers,  such  as  Clay  and  Medwin  [7],  etc.  Model 
tank  measurements  have  been  made  by  Kingsbury  [84],  Spindel  ]54] 
and  Zornig  [53],  etc.  A  relatively  simple  expression  that  fits  these 
measurements  with  a  fair  degree  of  realism  is 


2  2  2 

’•'(S.n.x)  -  exp{ - - ^ - -  }  .  cos  p(n-C  t)  (4.1.11) 

A  A  T.Z  P 

X  y  0 

This  represents  a  surface  corrugation  moving  in  the  y  direction  with 

a  velocity  C  .  The  parameters  A  ,  A  and  Tn  are  correlation  distances 
p  x  y  u 

in  the  x,  y  and  t  directions,  respectively.  This  expression  has  been 

used  in  scattering  studies  by  Clay  and  Medwin  [7],  McDonald  and 

Tuteur  [8].  A  shortcoming  of  equation  (4.1.11)  is  that  it  does  not 

take  into  account  the  surface  dispersion  which  is  an  obvious  and 

easily  observed  feature  of  typical  water  surfaces.  Dispersion  could 

be  modelled  by  using  a  stun  of  terms  of  the  form  given  in  equation 

(4.1.11)  with  different  values  of  p  and  C  ,  as  in  Tuteur,  McDonald 

P 

and  Tung  [40].  This  complication  has  not  been  added  to  our  analysis 
in  order  to  keep  the  result  simple. 


In  spite  of  the  relative  simplicity  of  equation  (4.1.11), 
the  expression  that  results  when  it  is  substituted  into  equation 
(4.1.6)  is  not  directly  integrable.  A  standard  procedure  that 
permits  integration  is  to  use  a  Taylor  series  expansion.  This 
can  be  applied  to  the  characteristic  function  to  give 


exp{-g2(l-'!'(C,n,T)]}  *  e-8  {l+g2f(‘)+  —  [g2'?(*)]2+ - } 


(4.1.12) 


(Dq 

where  g  =  —  Fa  is  the  Rayleigh  parameter  which  is  a  measure  of 
effective  surface  roughness  with  respect  to  acoustic  frequency .+ 


Since  |'P(£,n,x)|  <  1  ,  equation  (4.1.12)  converges  rapidly 
if  g  «  1;  hence  this  expansion  is  appropriate  if  the  surface 
roughness  is  small.  For  very  small  roughness  it  is  permissible 
to  use  only  a  first-order  approximation  of  equation  (4.1.12);  for 
surfaces  of  larger  roughness  (larger  g)  higher  order  terms  can  be 
used  as  well. 


For  a  very  rough  surface,  i.e.  g  »  1,  equation  (4.1.12) 
becomes  inconvenient  to  use.  In  this  case  we  observe  that  for 
very  large  g  the  expression  in  equation  (4.1.8)  has  a  significantly 
non-zero  value  only  when  l-’F(5,h,x)  is  almost  zero.  Therefore,  if 
we  expand  ¥(£;,n,x)  in  a  Taylor  series  about  5  ■  n  ■  x  *  0,  only 
the  first  few  terms  of  the  series  are  significant.  The  linear  or 
the  first  order  term  in  the  series  need  not  be  considered  since  it 


+  See  Chapter  II. 
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corresponds  to  a  non-differentiable  surface  [10],  [85].  Thus 
the  expansion  of  ’P(£,n,T)  correct  to  second  order  is 


«U,n,T)  «  1  - 


2  2  2 

(n-c  X)2 


(4.1.13) 


Using  this  second  order  expansion  of  ¥(£,11, x)  in  equation  (4.1.8)> 
we  get 

2  2  2  2 

exp{-g2[l-  »n ,t) ] >  *  exp{-  g2[  -4  +  JLr  +  -4  +  ^-{ii-C  t)2]} 

A  2  A  2  T  2  2  P 

x  y  0 

(4.1.14) 

The  other  statistical  parameters  that  need  to  be  specified 


are  the  slope  moments  in  J(s  )  .  The  orders  of  slope 


moments  in  J(?  )  range  from  1  to  4.  The  first  order  moment 

yl  y2 

is  the  average  surface  slope  which  is  zero.  The  second  order 

The 


2  *2  * 

moment  appears  in  the  form  of  C  ,  £  ,  and  C 


yly2 

first  two  are  the  mean  square  slopes.  We  assume  the  third  term 


to  ke  t*ie  form 


where 


C  ,C  =  E, 

yl  y2  2 


2  2 

exp{-  ~  -  ~  > 
A  t/ 

2  * 


(4.1.15) 


E2  * 


*2 


(4.1.16) 


clearly  £2  0  for  all  cases. 

The  third  order  slope  moment  appears  in  the  form 
2  *  *2 

C  C  +  C  5  ,  and  is  assumed  to  have  the  form 

yl  y2  yl  y2 


2 


*2 


(4.1.17) 
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where 


_2  r  3  -  -2  C  *3 

yl  y2 


(4.1.18) 


The  parameter  is  generally  positive,  and  the  appearance 

of  the  minus  sign  in  front  of  is  explained  in  Appendix  G.  Both 

the  second  and  the  third  order  slope  moments  are  functions  of  n 

2  *9 

and  x.  We  neglect  the  fourth  order  moment  t  £  ,  because 

yl  y2 

it  is  normally  much  smaller  and  it  does  not  contribute  any 
asymmetry  to  the  final  result.  This  follows  from  equation  (4.1,7). 
The  fourth  order  moment  comes  from  the  product  of  the  first  terms 
in  the  first  and  the  second  brackets,  which  is 


[  j  ^  F(l-£  2)][-j  ^  F(l-e  *2)]  =  (  —  F)2(1-2e  +  c  2  c  *2) 

cc  *y^  1  J  c  y2  c  2  vy^  y2  ' 


(4.1.19) 


Notice  that  the  mean  square  slope  term  is  replaced  by  its 

equivalent  e2<  It  will  be  shown  in  the  final  result  that  l-2e2 

represents  the  energy  left  in  the  coherent  component  because  of 

2  *2 

finite  surface  slopes.  X,  C  represents  a  higher  order 

yl  y2 

correction  term  to  this  coherent  response  only  and  therefore  can 


be  neglected. 
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4.2  Random  Surfaces  -  Small  Rayleigh  Parameter 


The  characteristic  function  for  a  Gaussian  surface  of 
Rayleigh  parameter  less  than  unity  was  given  in  equation  (4.1.12). 
For  the  small  Rayleigh  parameter  case,  only  the  first  few  terms 
are  important,  so  that  equation  (4.1.12)  is  rewritten  as 


QC^,^)  -  exp{j  F  ^(t)-j  F  ?2(t+t)> 


__c 

-g2  2  Ax2 

a  8  {B-g^[e 


cos  p  (n-C  t) ] } 
P 


(4.2.1) 


Details  of  the  evaluation  of  the  four-fold  integral  in 
equation  (4.1.6)  are  presented  in  Appendix  H;  this  integration  is 
basically  straightforward  but  lengthy.  The  result  can  be  put  into 
the  form: 


*(u0,t)  -  Kx[ I0  +  I ^  +  IXJ 


(4.2.2) 


where  definitions  for  IQ,  I1+  and  Ij_  are  given  in  equation  (H-9) , 
and  is  a  constant  coefficient. 

The  signal  spectrum  can  be  computed  by  substituting  equation 
(4.2.2)  into  equation  (4.1.4)  and  taking  the  Fourier  transform 
with  respect  to  the  variable  t. 


r(u,u0) 


00 

J"  r(t)r*(t+r)e  dx 
—00 


-  K, 


00 
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-j  (w-u>0)t 
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(4.2.3) 
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Details  of  this  computation  are  also  shown  in  Appendix  H.  The 
result  of  this  operation  is  contained  in  equation  (H-14) ,  repeated 
here  as  equation  (4.2.4): 


r(a),a>0)  -  K^SpQCw.WQj+Sp^Cu^Qi+Sp^Od.ioQ)]  (4.2.4) 


where 


SPQ 0*>,u)q)  -  f006(ui-u>0)+  ir*  E  Tmemfm0exp{-  “f-(o>-u>0) 2 }  (4.2.5a) 

m**2  ,3 


2  T0  2 

sp1+(<j,w0)  -  |-^{TofofexpI~  T~  (u_w6+  pCp)  J 


T  A  2(m-(o+pCp)2 

+  j  Vv  e  f  .  exp[ - 5 - 5 - ^ 

0^2,3  ^T02+Tm2)^  4(Tn2+T„2) 


TaT 
0  m 


k  .  ■  4. 

0  m 


(4.2.5b) 


The  signal  spectrum  is  seen  to  consist  of  three  parts; 
Sp0(u>,u>0),  Sp1+(a),u>0)  and  Spl^(u,a)Q).  Sp0(oj,co0)  represents  the 
coherent  response  of  the  scatter  channel.  Sp1+(w,10q)  and  Sp^_(o),oJq) 
are  respectively  the  upper  and  lower  shifted  Doppler  spectra.  The 
shape  of  both  spectra  is  approximately  Gaussian  with  center 
frequencies  at  uig+pCp  and  ui^-pCp.  The  two  center  frequencies  are 
the  transmitted  acoustic  frequency  Uq  shifted  up  and  down  by  pCp, 
the  frequency  of  the  surface  wave  correlation  function. 

The  width  of  SPq(u,u»q)  is  governed  by  Tj  and  T^,  which  are 
the  correlation  constants  (temporal)  of  the  second  and  third  order 


» 
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moments  of  surface  slopes.  The  widths  of  ^^(01,0)^)  and  S  ^  (u.Wq) 
are  controlled  by  the  combination  of  and  the  surface  wave- 

height  correlation  constant  Tq.  If  the  surface  correlation  functions 
are  in  the  form  of  an  undamped  sinusoid,  then  Tq,  T^  and  are 
infinite  and  the  signal  spectrum  consists  of  three  spectral  lines 
at  Wq,  Wq+pCp  and  ojQ-pCp .  This  is  essentially  the  same  result 
obtained  in  the  deterministic  case.  For  smaller  Tq,  Ty ,  and 
T^,  the  surface  wave  spectrum  is  itself  a  narrow  band  spectrum  and  the 
three  spectral  lines  spread  to  form  three  narrow  bands. 

There  would  have  been  more  than  just  one  Doppler  sideband  on 
each  side  of  the  carrier,  if  we  had  incorporated  more  terms  of  the 
expansion  of  the  characteristic  function  in  equation  (4.2.1).  The 
number  of  significant  sidebands  depends  on  the  magnitude  of  the 
Rayleigh  parameter  g. 

We  consider  next  the  relative  magnitudes  of  the  upper  and 
lower  sidebands.  As  in  the  deterministic  case,  the  results  depend 
on  wind  direction  and  we  therefore  consider  the  downwind  and  crosswind 
geometries  separately. 

Case  I.  Downwind  geometry 

The  computational  details  of  the  substitution  of  downwind 
geometry  parameters  in  the  signal  spectrum  are  presented  in  Appendix  H. 
The  result  is  expressed  in  equation  (H-24)  and  is  repeated  in  the 
following  equation 
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~2  tsD  (w,o>q)+Sd  (o»,o)q)+  Sd  (io,(Dq)]  (4.2.6) 


(roo+rior  0  1+  v  "1- 


where 

T  ^ 

22  l2  ,  .2 

-  p  0  e2  "  4 — <“-wo) 

Sjj  (uj,(o0)  ■  (l-2e2)<S(a)-w0)+(/rT  T2>*  - ^ — e 

0  g 


3  2 

poe.  -  t—  (u>-oi0) 

+(V?T  T,)  -  e 

3  g 


(4.2.7a) 
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*{(l-2€2)e  4  °  +  - — 

T02T32((o-a1(+PCp)2  8 
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±  —  •  ) 
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4(T02+T22) 


(4.2.7 .b) 


We  observe  Chat  the  coherent  response  Is  an  impulse  of 
(l-2e2)  2 

amplitude  - 5 —  e  6  at  frequency  oin.  If  the  surface  is 

<ro+ti> 

completely  flat,  i.e.  e2»e3=g*>0,  the  coherent  response  is 

the  only  term  in  rD(w,u>g)  and  it  has  the  form  (roo+rio^  ’ 

is  the  familiar  form  of  point  source  response.  When  the  surface 

2 

is  disturbed,  e2  and  g  are  nonzero.  The  attenuation  term  e  8 
results  from  the  surface  deformation  amplitude,  while  the  term 


202  represents  attenuation  caused  by  the  slope  terms.  To 

consider  the  effect  of  surface  waveheight  and  the  effect  of 

slopes  separately,  we  assume  the  surface  to  have  extremely 

small  deformation,  i.e.  extremely  small  g,  but  have  finite 

surface  slopes.  In  other  words,  the  surface  can  be  modelled 

as  a  group  of  plane  facets  located  on  the  flat  surface.  Each 

facet  is  random  and  has  its  own  slope  and  the  correlation  of 

the  slopes  is  described  by  equations  (4.1.15)  and  (4.1.17). 

The  spectra  S  (w,mn)  can  be  neglected  for  extremely  small  g 
_  2  1±  ° 

and  e  8  is  approximated  by  1.  The  spectrum  in  equation  (4.2.6) 
can  then  be  rewritten  as 


rD(“-V  ■  (rJtiy(1-2e2><("-"oH(VT  T2> 


2  2  2  ,  .2 
p  a  e2  -  — (w-Up) 

- ^ —  e 


l00  ‘10 


poe3  - 
+  (*T  T3)  e 


3  <  \2 


(4.2.8) 


It  is  seen  that  the  energy  which  is  removed  from  the  coherent 
response  is  distributed  in  the  frequencies  surrounding  i.e. 
the  second  and  the  third  terms  in  the  brackets  in  equation  (4.2.8). 
As  discussed  in  the  previous  section,  the  bandwidth  is  determined 
by  Tj  and  Tj,  and  is  not  a  function  of  slope  moments.  Varying  the 
mean  square  slopes  ^2  only  changes  the  absolute  magnitude  of  each 
spectral  component  proportionally.  Figure  4.2  contains  sketches  of 


rD(u>,w0)  of  equation  (4.2.8)  for  both  small  and  large  values  of  £2' 
Notice  that  the  bandwidth  stays  the  same  in  both  cases. 
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Figure  4.2  Power  spectrum  -  small  and  large  e2 
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The  effect  of  somewhat  larger  surface  deformation  on 

the  frequency  spreading  function  can  be  examined  by  including 

the  two  spectra  S  and  S_^  (to,u>n)  into  our  considera- 

1+  u  °1-  u 

tion.  Each  of  the  spectra  is  the  sum  of  three  Gaussian  functions 

as  indicated  by  equation  (4.2.7b).  Considering  the  bandwidths 

of  these  three  Gaussian  functions  in  more  detail,  we  see  that 

they  are  respectively  2 (it  Tq)  2:r  ^(Tq~^+  T2  2)^  and 
-1  -2  -2  k 

2u  (Tq  +T^  )  Hz.  Hence  the  spread  in  the  first  Gaussian  is 

due  entirely  to  the  surface  waveheight  fluctuations,  while  the 
spreads  in  the  second  and  the  third  Gaussian  functions  are  due 
to  the  combined  effects  of  surface  f luctuations  and  surface 
slopes.  The  combination  is  such  that  the  bandwidth  is  controlled 
by  the  process  which  has  the  smaller  correlation  time.  To  our 
knowledge,  no  measurements  of  slope  correlation  times  have  ever 
been  made.  To  be  consistent  with  the  assumption  that  slopes  and 
heights  are  independent  because  the  slopes  depend  largely  on 
small  ripples  one  might  argue  that  T2  T^  <<  Tq  .  However, 
another  point  of  view  is  to  assume  T2 «  T^  k  Tq,  simply  on  the 
basis  that  all  three  correlation  times  result  from  the  same 
surface  wave  motion.  We  have  used  the  second  approach  in  equation 
(4.2.7b)  to  compute  Sp  (uj,u)q)  and  Sp  (to, Wq)  .  All  the  other 
parameter  values  needed  for  this  calculation  are  taken  from  data 
obtained  in  a  model  tank.^  These  values  are  shown  in  Figure  4.3, 
which  contains  a  plot  of  the  two  spectra.  The  Rayleigh  parameter 

^  The  magnitudes  of  c2  and  used  in  the  plot  may  be  larger  than 
the  practical  value  in  order  to  emphasize  the  asymmetry. 
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g  is  0.3.  We  observe  that  the  two  sideband  spectra  have 


completely  different  amplitudes.  The  amplitude  differences 

pcre. 

are  caused  by  the  two  double-signed  terms,  i.e.  + -  and 

2  2  6 
4o  R  dp 

+  ■-  9 -  .  Note  that  the  inequality  of  the  spectral 

g  L  (m  ) 

heights  does  not  depend  on  the  correlation  times  and  is  therefore 
unaffected  by  the  assumptions  made  about  the  relative  values  of 
Tq,  and  T^.  As  in  the  deterministic  case  considered  in 
Chapter  III,  the  two  double-signed  terms  can  be  related  to 
two  different  mechanisms,  i.e.  the  surface  slope  asymmetry 
and  the  scattering  geometry  asymmetry.  The  effect  of  surface 


slope  asymmetry  shows  in  the  term 


poe, 

g 


,  while  the  effect  of 
2„2 . 


4o  R  dp 

scattering  geometric  asymmetry  is  in  — -  . 

gVoo0) 

We  first  consider  only  the  surface  slope  asymmetry  by 
setting  d  «  0.  From  equations  (H-13)  and  (H-22),  this  implies 
either 


T  m  T 
00  10 


(4.2.9a) 


or 


* 0  *  90 


(4.2.9b) 


Since  e^>  0  for  wind-driven  surface,  we  clearly  see  that  the 
upper  sideband  S  (w,u>n)  is  larger  than  the  lower  sideband 
Sp  (u,Wq).  If  the  direction  of  surface  wave  propagation 
reverses,  both  Cp  and  change  sign  so  that  the  sideband 
asymmetry  stays  the  same.  The  degree  of  this  sideband  asymmetry 
is  a  function  of  acoustic  frequency.  If  we  increase  the  acoustic 
frequency,  the  Rayleigh  parameter  g  becomes  larger,  and  the  term 
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p<7£3 

— — —  becomes  smaller.  Therefore,  we  would  expect  to  find  less 
sideband  asymmetry  for  higher  u>q.  This  is  illustrated  in  the  two 
sketches  of  Figure  4.4.  The  acoustic  frequency  used  to  generate 
the  lower  sketch  is  twice  the  value  used  in  the  upper  sketch. 

We  now  consider  the  sideband  asymmetry  caused  by  asymmetric 
scattering  geometry.  In  order  to  do  this,  we  let  ■  0  to 
elminate  the  effect  of  surface  slopes.  The  form  of  the  exponential 
term  can  be  rewritten  by  replacing  d  with  its  equivalent  expression 
in  equation  (H-22) .  We  get 


exp{  + 


4o2R2dp 
g2L2 (w0) 


P 

exp  {  +  [  - 

“ o/c 


rio~roo 

r10+r00 


2  cos  i/>n 

- - - ±  ]}  (4.2.10) 

sin  i|»0 


This  is  similar  to  equation  (3.2.18)  for  very  large  L(w0).  The 
upper  and  the  lower  sidebands  are  asymmetric  when  rQQ  4  and 
i  90°.  The  discussion  of  the  sideband  asymmetry  is  the  same 
as  in  Chapter  III.  Figure  4.3  and  4.6  contain  sketches  of  the 
asymmetric  sidebands  for  Tqq  >  r^Q  and  Tqq  <  r^Q  ,  respectively. 
The  corresponding  downwind  geometry  is  shown  above  each  sketch. 


Case  II .  Crosswind 

The  formulation  of  the  signal  spectrum  in  a  crosswind  geometry 
is  presented  in  Appendix  H.  The  resulting  spectrum  is  written  as 


r(u,u0)  -  - - - r  [S  (u,o>0)+  S  (w,o>q)+  Sc  (w.Uq)]  (4.2.11) 

(roo+rio)  0  1+ 
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DIRECTION 


Z 


Asymmetric  sidebands  for  small  Rayleigh  parameter 
e3  =  0,  and  rQ0  >  r1Q 


4-24 


<r-"  ' 


where 


2  2  2  2 
p  o  e  -  T—  (w-w0) 

Sc  (to,uQ)  -  (l-2e2)6(a)-ui0)+(V:!r  T2>  — ■£%-=■  e 


+(</"  T3)  — —  e 


3  2 
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(4.2.12b) 


We  observe  that  equations  (4.2.12)  are  quite  similar  to 


equations  (4.2.7)  except  that  the  geometry-asymmetry  factors, 

2  2 

i.e.  exp{+  4yy--'P—  }  in  equation  (4.2.7b),  have  disappeared. 

gV(u>0) 

Therefore,  the  siaeband  ratio,  i.e,  Sc^+(w,u)q)/Scj_(u),<i)q)  , 
will  remain  unchanged  for  all  values  of  r^  and  r^Q.  The  effect 
of  asymmetric  slope  distribution  is,  however,  still  present  in 
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equation  (4.2.12).  Figure  4.7  shows  a  crosswind  spectrum  obtained 
by  Zornig  [3]  in  a  model  tank  experiment. 


Figure  4.7  Experimentally  measured  spectrum-crosswind,  wind  speed  at  5.4  m/sec 
f n  =  524  KHz,  r,,n  —  r._  **  173  cm,  =  17 


4.3  Random  surfaces  -  large  Rayleigh  parameter 


When  the  surfaces  become  rougher,  i.e.  g  =  — - —  ;»  X,  either 
by  increasing  the  r.m.s.  surface  waveheight  a  or  by  increasing  the 
acoustic  frequency  u^;  the  surface  height  characteristic  function  is 
conveniently  approximated  by  equation  (4.1.14)  repeated  here  as 
equation  (4.3.1): 


Q(C1.?2^  *  exp{-g2[  -£j  +  +  £-  (n-Cpx)2]} 


(4.3.1) 


The  received  signal  spectrum  can  be  computed  as  in  Sec.  4.2, 
with  equation  (4.3.1)  being  used  to  represent  the  characteristic 
function  Appendix  I  has  the  details  of  this  computation. 

The  resulting  spectrum  is  expressed  in  equation  (1-4) ,  repeated 
here  as  equation  (4.3.2): 
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(4.3.2) 
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where  K2>  f1Q,  ^3*  ^  and  ^2  are  defined  in  Appendix  I. 

Notice  that  the  signal  spectrum  r(u>,u>Q)  is  the  superposition 

of  three  functions.  Each  of  the  three  functions  contains  a  Gaussian 

term.  The  centers  of  the  first,  the  second,  and  the  third  Gaussian 

are  respectively  +  q2(0),  +  q2(  )»  and  “q  +  q3^  ' 

A 2 

From  equation  (I-4b) ,  we  get 
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(4.3.3b) 
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(4.3.3c) 


where  N  is  a  complex  quantity  defined  in  equation  (H-4) .  Equation 
(4.3.2)  is  rather  unwieldly  as  it  stands  and  must  be  simplified 
considerably  to  permit  easy  interpretation.  We  consider  first  the 
equation  for  q2(").  From  equation  (4.3.3)  we  see  that  their  dependence 
on  A2  and  A^  is  relatively  slight,  especially  if  we  make  the  assumption 
that  A2>  A^  and  A^  are  all  of  the  same  order  of  magnitude  (this 
assumption  is  based  on  considerations  similar  to  those  used  in  the 
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previous  section  to  argue  that  T^~  « 

2 

dealing  with  rough  surfaces  here,  g  » 
Thus  we  have  that 


T^).  Note  that  since  we  are 

2  2  11 

1,  and  therefore  g  /A  >>  — r-, — 7 

y  a  1  a  *■ 


q 


2 


q2(0) 


(4.3.4) 


Furthermore,  if  the  surface  slope  correlation  times  and  are  of  the 
same  order  of  magnitude  as  Tq,  then  from  equation  (I-4a) ,  we  have 
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75  >  - 

A„ 


(4.3.5) 


By  the  use  of  equations  (4.3.4)  and  (4.3.5)  the  three  Gaussian 
functions  in  equation  (4.3.2)  are  made  to  be  identical  and  r(u>,w0)  can 
be  rewritten  as 


r(w,wQ) 
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TT 

q^o.o) 


h 

)  • 


[w-w  -q2(0)]2 

-p{-  -  4qikoy 


-j[o)-u)Q-q  (0)]  -q  (0)] 

{f10+  e2f12(  2q1(0>0)  )  +  e3f13(  2^(o7o)  ) 

(4.3.6) 

Except  for  the  terms  involving  the  slope  moments  e 2  an d  e^*  the  signal 

spectrum  r(w,Wg)  is  seen  to  be  a  Gaussian  function  with  center 

-1  U 

frequency  at  +  c^(0)  and  bandwidth  ffi  ■  2ir  [q^(0,0)]  .  The 
functions  f^Q,  f^,  f^»  q^  and  q^  all  depend  on  the  scattering 
geometry,  as  cai  be  seen  in  Appendix  I.  Thus  the  spectrum  r(u),u)g) 
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depends  on  geometry  and  on  slope  moments.  As  in  the  previous  sections 
we  consider  in  detail  only  the  down  (up)  wind  and  crosswind  geometries. 


» 


» 


» 


I 


Case  1.  Downwind 

Detailed  calculations  for  the  downwind  case  are  given  in 
Appendix  I,  equations  (1-5)  to  (1-9).  Equation  (1-9)  gives  the 
expression  for  the  spectrum,  repeated  here  as  equation  (4.3.7) 
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We  see  that  r(u),u)g)  is  generally  not  symmetric  with  respect  to  u)Q 
unless  ^2(0)  “  0  and  ■  0.  These  are  the  effects  of  the  two  mechanisms 
mentioned  previously  as  we  recognize  that  represents  the  slope  asymmetry 
and  q2(0)  is  the  scattering-geometry  asymmetry.  The  expression  of  q2(0) 
can  be  obtained  from  equation  (4.3.3a)  and  the  expression  of  d  in  equation 
(H-22) 
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We  see  that  q2(0)  is  non-zero  if  J*  r^g  and  90°;  and 

q2(0)  is  positive  (negative)  when  r^g  <  rig(rQQ  >  rl0  ^  •  Therefore, 
considering  no  slope  asymmetry,  i.e.  “  0,  we  conclude  that  if  the 
wind  blows  from  source  to  receiver,  and  if  r^  <  r1()(rQQ  >  r1Q)  »  the 
upper  (lower)  sideband  is  larger  than  the  lower  (upper)  sideband 
(see  Figure  4.8).  This  conclusion  is  the  same  as  was  reached  in  the 
deterministic  case  and  in  the  slightly  rough  (random)  surface  case. 

We  examine  next  the  sideband  asymmetry  caused  by  surface  slope 
asymmetry.  As  usual,  we  eliminate  the  effect  of  the  scattering  gepjjetry 
factor  by  assuming  r^g  =  r10*  T^en  equation  (4.3.7)  becomes 

2 

i  J —  (w-i O 

r(u»,a)0)  »  - r  •  - —  *  expf - r-  ) 

(r00+r10)  gpCp/2  (2irfB) 

j  ,2  OE- 

•{(l-2e_)+2o  e„[(  — -  +  )+  — - — t-  .]+ -  (w-u>  )}  (4.3.9) 

\  2  g2A/  gCp  ° 


We  see  that  the  center  of  the  Gaussian  function  is  at  Ug.  The 

sideband  asymmetry  is  introduced  by  the  last  term  in  the  brackets,  viz. 
°e3 

(u-Wq) ,  which  is  positive  and  therefore  Increases  the  upper  sideband 
when  u>Wg;  and  is  negative  and  decreases  the  lower  sideband  when  o)<o)q. 
Therefore,  the  upper  sideband  is  always  larger  than  the  lower  sideband. 


Figure  4.9  contains  a  plot  of  equation  (4.3.14)  with  fg  **  256  KHz, 

^g  ■  30°,  rQQ  *  r^g  ■  254  cm.  A  slight  asymmetry  can  be  seen  by  comparing 
the  corresponding  spectral  amplitudes  on  both  the  upper  and  lower  sidebands. 
Much  greater  asymmetry  is  produced  by  differences  in  source  and  receiver 
depths.  This  is  illustrated  in  Figures  4.10  to  4.12,  which  are  plotted 
for  roo/rio  “  °*1*  10,  and  0.5  respectively. 
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Figure  4.8  Gaussian  spectrum-Downwind 


o 


4-35 


Figure  4.9  Theoretically  computed  power  spectrum  (1)  -  rough  surface,  upwind 
Parameters  used  are  rQ0  =  r  =  254  cm,  =  30°,  o  =  0.2  cm, 
fn  =  256  KHz,  Cp  =  30  cm/sec,  T_  =  0.1  sec,  e9  =  0.3,  =  0.2. 
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The  amount  of  frequency  spreading  is  determined  by  the  bandwidth 
From  equation  (I-4a) ,  with  a  ■  b  ■  0,  we  find  that  this  is  given 
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(4.3.10) 


where  the  Rayleigh  parameter  g  has  been  replaced  by  g  = 


uj0Fo  (2irf0)Fo 


We  see  that  the  bandwidth  fg  is  a  monotonic  function  of  the 

acoustic  frequency  f^  (in  Hz).  However,  it  is  not  a  linear  function  of 

f  ,  since  the  beamwidth  L(a)_)<*  is  frequency  dependent.  For  low 

0  u  £o 

acoustic  frequency,  L(ojq)  is  large,'  and  we  can  assume  that 
L2(m0)  1  2 

— rr~  »  — r  +  ±—  .  Hence  for  small  £n 

8  a¥  A  1  2  U 


(4.3.11) 


As  the  acoustic  frequency  increases,  the  beamwidth  L(u«)  becomes  smaller, 

l/(cO  1 


L<“0> 

and  at  some  high  frequency,  we  have  — 5—= — 

8a  V 

fn,  e*luation  (4.3.10)  can  be  rewritten  as 


Hence,  for  large 


£  s  f  /is  r  fV  1  ) 
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Figure  4.13  contains  a  plot  of  fg/fg  as  function  of  frequency  fg 
according  to  equation  (4.3.10).  Notice  the  asymptotic  behavior  at  both 


low  and  high  acoustic  frequencies. 


Case  II.  Crosswind 


The  received  signal  spectrum  for  crosswind  is  derived  in  Appendix 
I,  equation  (1-10),  repeated  here  as  equation  (4.3.13) 
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{(1-2c2)«»2c2  [  <  A  +  f  >+  XT ]+  °T  (4-3-13) 

A  2  g  A.  gCp  7 


There  is  only  one  factor  in  equation  (4.3.16)  that  can  cause  side- 
ae3 

band  asymmetry.  It  is  (o>-U)q)  ,  the  surface  slope  asymmetry.  As 

we  discussed  in  the  section  dealing  with  the  downwind  geometry,  this 
term  will  make  the  upper  sideband  larger  than  the  lower  sideband  in  a  wind 
driven  surface.  Varying  the  scattering  geometries,  i.e.  changing  rQ0,  r1Q, 
has  no  effect  on  the  sideband  asymmetries. 


The  bandwidth  for  crosswind  is  obtained  by  substituting  crosswind 


parameter  values  into  equation  (I-4a).  This  gives 


l/(u>0)  1 

r  1  .  pV  a/ 

2  +  2 

Lt0  2  L  <«n>  1 

2  J  f 

P 

(4.3.14) 


2  2  2 
8o  R  sin 


Very  low  values  of  f^  might  be  Inconsistent  with  the  assumption  of  a  very 
rough  surface;  i.e.  g  >>  1.  However  the  inequality  leading  to  equation 
4.3.11  depends  on  beam  width,  which  is  partly  controlled  by  transducer  size. 
Hence  equation  4.3.11  may  hold  with  g  >>  1. 
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The  main  difference  between  equations  (A. 3. 10)  and  (A.3.1A)  is  that 


the  term 


l2(uq) 

8o2R2 


in 


L2 


in¬ 


equation  (A. 3. 10)  is  replaced  by  — 7  s  ~ 


in  equation  (A.3.1A).  Since 


l2(V 


8a  R  Sin  i|)q 
is  always  larger  than 


L2((o0)  8a2R2Sin2^0 

,  the  crosswind  bandwidth  in  equation  (A.3.1A)  is  larger  than 


8c2R2 


the  downwind  bandwidth  in  equation  (A. 3. 10).  In  other  respects  the 
dependence  of  the  bandwidth  on  fg  is  similar  to  that  already  discussed 
for  the  downwind  case.  Figure  A.1A  shows  plots  of  bandwidths  for  both 
downwind  and  crosswind  at  ij/g  =  17°  and  30° .  Also  shown  in  this  figure 
are  experimental  points  of  bandwidth  measurements  made  in  a  model  tank 
by  J.G.  Zornig  [3]  .  The  general  behavior  of  the  theoretical  and  experi¬ 
mental  results  are  very  similar,  and  the  correspondence  between  the 
results  for  ijjg  =  30°  are  well  within  experimental  error.  For  iJ>q  =  17° 
there  may  already  be  some  shadowing,  and  this  might  explain  the  con¬ 
sistently  smaller  bandwidth  observed  in  the  experiment.  This  issue  has 
not  yet  been  investigated. 


In  Figures  A. 15  to  A. 17  are  plots  of  the  crosswind  signal  spectrum 
according  to  equation  (A. 3. 13);  in  Figure  A. 15,  fg  =■  256  KHz  and  \|>g  =  17°; 
in  Figure  A.  16,  fg  =  256  KHz,  i)jg  =  30°,  and  in  Figure  A.  17,  fg=  52A  KHz 
and  ipg  =  30°.  Note  that  the  sideband  asymmetry  becomes  less  significant 
with  increase  in  fg.  Also  note  that  the  bandwidth  increases  both  with 
i|ig  and  with  fg. 

In  summary,  we  have  examined  the  frequency  spreading  function  of 
the  acoustic  signal  scattered  from  a  random  rough  surfaces  for  both  the 
slightly  rough  case  and  very  rough  case.  We  have  found  that  the  amount 
of  frequency  spreading  (bandwidth)  depends  on  both  the  surface  roughness 
and  surface  slopes.  The  following  illustrations  summarize  the  results  of 


section  A. 2  and  A. 3 


The  incident  acoustic  signal  is 
sinusoidal  with  frequency  u q. 


1  w 

“0 

The  received  signal  spectrum  re¬ 
flected  from  an  absolutely  flat 
surface.  The  spectral  amplitude 
is  attenuated. 

U) 


Sinusoidal  surface  of  frequency  ft. 
The  scattered  signal  spectrum  is 
a  sequence  of  spectral  lines  of 
Interval  ft. 


I  1  1  1  1  I  1.4 

■4fi  «0-2fi  (i)q  0JQ+2R  Wq+AQ 

-*1 


0) 
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Surface  with  very  small  roughness 
finite  surface  slopes  with  large 
temporal  correlations. 


J 


,■4*-  .«* 


Same  as  (A) ,  except  the  slope 
has  smaller  temporal  correlation 
(become  more  uncorrelated) . 


From  (A) ,  the  surface  roughness 
Increases,  and  the  surface  wave 
Is  a  narrow  band  process  with  a 
center  frequency  of  pCp. 


The  surface  roughness  Increases 
more,  and  the  bandwidth  of  the 
surface  wave  also  Increases. 


A  continuation  of  (7) . 


For  very  large  surface  roughness, 
the  signal  spectra  merge  Into  one 


Gaussian 


The  amplitudes  of  the  frequency  spreading  function  are  in  general 
not  symmetric  with  respect  to  the  transmitted  frequency.  There  are  two 
different  mechanisms  which  are  responsible  for  this  amplitude  difference. 

One  mechanism  is  the  surface  slope  asymmetry,  i.e.  is  positive,  which 
always  results  in  the  upper  sideband  being  larger  than  the  lower  sideband. 

This  effect  does  not  depend  on  the  wind  direction.  The  other  mechanism  is 
the  scattering  geometry  asymmetry  caused  by  unequal  depths  of  source  and 
receiver.  This  effect  is  a  function  of  the  wind  direction  and  depends 
on  the  relative  depths  of  source  and  receiver.  Specifically,  if  the 
source  is  located  at  a  smaller  depth  and  if  the  wind  blows  from  source 
to  receiver,  then  the  upper  sideband  is  larger  than  the  lower  sideband. 

By  either  reversing  the  wind  direction  or  maintaining  the  wind  direction 
but  placing  the  receiver  at  a  smaller  depth  than  the  source,  the  lower 
sideband  will  have  a  larger  amplitude  than  the  upper  sideband.  Since 
this  mechanism  is  wind-direction  dependent,  it  has  no  effect  on  the  side¬ 
band  amplitude  asymmetries  in  the  crosswind  condition.  The  effect  of  the 
second  mechanism  on  the  sideband  asymmetry  is  in  general  larger  than  that 
of  the  asymmetric  surface  slopes.  The  effect  of  both  mechanisms  decreases 
as  the  acoustic  frequency  increases.  In  other  words,  the  frequency  spreading 
function  is  more  symmetric  at  higher  acoustic  frequency  than  at  lower 
frequency. 

The  width  of  the  frequency  function  is  controlled  by  several  parameters; 
it  increases  with  the  grazing  angle,  the  surface  roughness,  and  the  acoustic 
frequency.  It  also  depends  on  the  wind  direction,  the  width  of  the  spreading 
increases  as  the  geometry  changes  from  downwind  to  crosswind. 
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Chapter  V  Experimental  Data  Comparison  and  Discussion 
5.0  Introduction 

A  series  of  surface-scatter  experiments  has  been  performed  in 
a  model  tank  by  J.G.  Zornig  [3]  .  A  windblown  surface  is  generated  in 
this  tank  by  a  fan  blowing  air  through  a  wind  tunnel  suspended  above 
the  surface.  A  detailed  description  of  the  experimental  equipment 
and  setup  can  be  found  in  [79].  The  experimental  technique  involves 
the  transmission  of  narrow  pulses  that  are  preemphasized  so  as  to  produce 
an  approximately  flat  received  spectrum  over  the  frequency  range  of  0.2 
to  1.2  MHz.  The  received  pulses  are  acquired  by  a  digital  computer  and 
stored  on  magnetic  tape.  Data  processing  to  obtain  the  frequency  spreading 
function  can  then  be  performed  on  the  stored  data.  The  data  processing 
techniques  are  described  in  detail  in  Appendix  J. 

5.1  Doppler  Shift 

Examples  of  measured  and  computed  spectra  are  shown  in  Figures 

5.1  to  5.4.  Each  figure  contains  four  spectra  corresponding  to  four 
transmitted  frequencies,  f q  :  256  KHz,  524  KHz,  792  KHz,  and  1.06  MHz. 

The  zero  frequency  of  all  plots  is  referred  to  the  transmitted  frequency. 
Hence  the  positive  (negative)  frequencies  correspond  to  the  upper  (lower) 
sideband.  Figures  5.1  and  5.3  are  for  upwind  transmission,  while  Figures 

5.2  and  5.4  are  for  crosswind.  The  grazing  angle  in  Figures  5.1  and  5.2  is 
30°  and  in  Figures  5.3  and  5.4  it  is  17°.  In  all  plots  the  Jepths  of 
source  and  receiver  are  the  same.  It  is  clear  from  these  figures  that 
there  is  a  significant  discrepancy  between  theory  and  experimental  results. 
According  to  the  theory,  the  asymmetry  should  decrease  with  f^. 


This 


follows  from  equations  (4.3.7)  and  (4.3.13).  The  term  that  causes 
asymmetry  Is  the  term,  and  Its  coefficient  is  inversely  proportional 
to  acoustic  frequency  Uq.  However,  the  discrepancy  between  theory  and 
experiment  is  most  significant  at  high  acoustic  frequencies.  It 
appears  in  all  of  these  figures  that  there  is  a  doppler  shift  in 
addition  to  the  theoretically  predicted  side-band  asymmetry. 

A  possible  explanation  for  this  doppler  shift  can  be  found  in 
small  azimuthal  misalignments  of  the  transducers  used  in  the  experiment. 
Even  though  the  transducer  mountings  were  carefully  positioned  to  be 
exactly  opposite  each  other  with  respect  to  the  specular  point,  the 
possibility  remains  that  small  imperfections  in  the  transducers  them¬ 
selves  could  cause  the  center  lines  of  their  beam  patterns  to  be  a  few 
degrees  off  the  boresight  direction.  Such  a  divergence  is,  in  fact, 
indicated  in  [53]. 

To  test  this  conjecture  several  experiments  were  performed. 

The  results  of  these  experiments  are  shown  in  Figures  5.5  to  5.7. 

Figure  5.6  was  obtained  from  a  crosswind  geometry  with  <(>T  ■  270°  and 
$R  ■  90°,  while  Figure  5.5  ($T  *>  270°,  <j>R  *  85°) and  Figure  5.7  ($T  =  270°, 
$R  =  95°)  were  obtained  in  a  roughly  crosswind  geometry  but  with  receiver 
being  intentionally  misaligned  by  5°.  The  top  view  of  the  relative 
source  and  receiver  positions  in  these  experiments  is  shown  in  Figure 
5.8.  In  this  figure,  the  source-receiver  pairs  for  Figures  5.5,  5.6  and 
5.7  are  respectively  (S,  R^) ,  (S,  R2)  and  (S,  R^) .  Since  in  Figure  5.6, 
the  source-receiver  pair  (S,  was  supposedly  exactly  crosswind,  we 
would  expect  a  roughly  symmetrical  spectrum  for  high  acoustic  frequency. 
However,  we  find  that  the  most  symmetrical  spectrum  among  these  figures 
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Power  spectrum  -  experimentally  measured  at  $  =  270 


is  contained  in  Figure  5.5,  i.e.  for  the  source-receiver  pair  (S,  R^) . 

In  Figure  5.6,  we  see  a  fair  amount  of  frequency  shift,  approximately 
12  Hz.,  in  the  plot  for  f^  =  1.06  MHz.  There  is  a  even  larger  frequency 
shift  in  Figure  5.7  almost  45  Hz  toward  the  positive  side  in  the  plot 
for  fg  =  1.06  MHz. 

If  we  assume  that  the  center  lines  of  the  source  and  receiver 
beam  patterns  diverge  from  their  respective  bore-sight  direction  by 
2.5°,  then  it  is  conceivable  that  the  intentional  misalignment  in 
Figure  5.5  could  result  in  a  roughly  crosswind  geometry,  with  the 
specular  point  at  O'  as  shown  in  Figure  5.9.  The  other  two  gemetries 
would  then  represent  angled  paths,  with  both  source  and  receiver 
pointed  slightly  upwind.  Since  both  beam  patterns  would  then  be  larger 
in  the  upwind  direction  than  in  the  downwind  direction,  we  would  expect 
more  energy  returned  with  upshifted  frequency.  The  frequency  displace¬ 
ment  should  be  larger  for  the  receiver  at  R^  than  for  the  one  at  R2» 

This  is  exactly  what  is  observed  in  Figures  5.5  to  5.7 

5.2  Theoretical  Analysis  of  Angled  Paths 

A  theoretical  treatment  of  surface  scatter  over  a  slightly  angled 
path  is  easily  obtained  by  modifying  the  beam  pattern  function  used 
in  the  theory.  Assume  that  the  center  of  the  beam  pattern  function  on 
the  water  surface  is  off  the  specular  point  by  a  distance  Ax  in  the  x 
direction  and  Ay  in  the  y  direction.  Then  the  form  of  the  modified 
beam  function  can  be  derived  from  equation  (2.3.13): 
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where  L  (u>.) ,  L  (w.)  and  L  (w_)  are  defined  in  equation  (2.3.14) 
xx  O’yyO  xy  0 

Rearranging  the  terms  in  the  brackets  results  in 
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The  first  three  terms  inside  the  brackets  represent  the  original 
beam  function.  The  last  three  terms  represent  an  additional  attenuation 
of  the  received  acoustic  power  due  to  the  off-specular  beam  function. 
They  don't  have  any  effect  on  the  sideband  asymmetry  or  the  doppler 
shift  and  therefore  can  be  neglected.  The  two  remaining  terms  represent 
the  significant  modifications  of  the  theory.  Substitution  of  the  new 
beam  pattern  function  of  equation  (5.2.2)  in  equation  (F-9)  affects 
mainly  the  coefficients  D  and  E.  If  we  represent  the  new  values  as 
and  E^,  then  we  have 
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The  shift  of  the  center  of  the  Gaussian  spectrum  was  defined  as 
q2(0)  in  equation  (4.3.3a).  The  parameter  d  appearing  in  this  equation 
is  defined  in  equation  (H-13)  and  is  seen  to  depend  on  E  in  equation 
(H-22).  If  we  use  instead  of  E  in  this  definition,  q3(0)can  be 
rewritten  as  follows: 
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We  see  that  f^  is  the  additional  doppler  shift  (in  Hz)  introduced  by 

the  misalignment.  As  we  assumed  in  the  previous  section,  the  center 

of  the  beam  pattern  function  is  at  0’  as  in  Figure  5.9.  Thus  the  beam 

center  displacements  are  A  ■  0  and  A  ■  -  1 00* | ,  where  1 00 T |  indicates 

x  y 

the  distance  between  0  and  O'.  Hence  the  frequency  displacement 


1 0011  a  positive  quantity,  i.e.  the  Gaussian  spectrum 


M  -  2,  s 

*  Lyy 

is  displaced  toward  positive  frequencies.  This  is  what  was  observed  in 


Figures  5.5  to  5.7.  In  fact,  the  frequency  displacements  shown  in  all 
experiments  shown  in  Figures  5.1  to  5.4  can  be  explained  in  this  way. 


The  dependence  of  on  the  acoustic  frequency  fg  can  be  computed 
by  using  equations  (2.3.14)  and  (5.2.6a)  in  equation  (5.2.6b).  An 
expression  for  B,  to  be  used  in  equation  (5.2.6a)  can  be  found  in 
Appendix  F,  equation  (F-10) .  When  this  is  done  we  find  that 

2  2r 
g  P  Cp 

fM  =  j“jl  2  2~2  '  7"  ‘  ^OO^IO’V^T'V 

|N|  g  P  g 

N+N*  A  2  2  (5.2.7) 

y 

where  h(rpQ,r^,^T>^^,4i^,4>R)  is  a  function  of  geometric  parameters  only. 

We  show  in  equation  (H-20)  that  for  downwind  for  crosswind  the  term 
2  u>02l2(u)0) 

|N]  /(N+N*)  can  be  approximated  by  - t—t —  if  the  beam  pattern 

2c  K 

illuminates  several  Fresnel  zones.  Since  the  beam  pattern  width 

2  2 

is  proportional  to  1/u)q,  (Oq  L  (wq)  is  independent  of  frequency,  and 

2  2 
therefore  |n|  /(N+N*)  is  a  constant.  On  the  other  hand,  g  is  pro- 

2  2 
portional  to  Uq  .  For  frequencies  high  enough  to  make  |n|  /(N+N*) 

negligible  in  the  denominator  of  equation  (5.2.7),  fM  is  proportional 

to  <Uq.  Figure  5.10  contains  a  plot  of  the  measured  frequency  displacements 

for  the  experimental  results  shown  in  Figure  5.5  to  5.7.  The  linear 

dependence  between  f„  and  f.  is  clearly  shown. 

M  U 

Figures  5.11  to  5.14  show  the  same  experimental  spectra  as  those 
shown  in  Figures  5.1  to  5.4  and  a  set  of  theoretical  spectra  shifted  by 
values  of  f^  chosen  to  provide  the  best  possible  match  with  f^  constrained 
to  be  linearly  dependent  on  f^.  These  values  are  shown  in  each  figure. 

It  is  seen  that  the  correspondence  between  the  theoretical  and  experi¬ 
mental  results  is  very  good. 


Frequency  displacements  for  spectra  shown  in  Figures  5.5  to  5.7. 


Shifted  power  spectrum  of  Figure  5.1.  The  frequency  displacements  for  fp  =  256  KHz,  524  KHz 
792  KHz,  and  1.06  MHz  are  respectively  3  Hz,  6.1  Hz,  9.1  Hz,  and  12  Hz. 


Shifted  power  spectrum  of  Figure  5.2.  The  frequency  displacements  for  f~  *  256  KHz 
524  KHz,  792  KHz,  and  1.06  MHz  are  respectively  2  Hz,  4  Hz,  6  Hz,  and  8  Hz. 


5 . 3  Bandwidth 


According  to  equations  (4.3.10)  and  (4.3.14),  the  bandwidth  should 
increase  with  acoustic  frequency  f^,  with  grazing  angle  \|/q,  with  surface 
roughness,  and  with  size  of  beam  pattern.  The  increase  in  bandwidth  with 
f^  appears  in  all  of  the  figures;  the  effect  of  grazing  angle  increase 
can  be  seen  in  a  comparison  of  Figures  5.11  and  5.13,  also  in  Figures 
5.12  and  5.14;  the  effect  of  beam  pattern  increases  between  Figure  5,12 
(small  band  width)  and  Figure  5.6  (large  beam  width)  is  caused  by  increasing 
r^Q  and  r^Q.  An  additional  qualitative  check  of  the  theory  is  in  the 
bandwidth  comparison  between  upwind  and  crosswind;  the  theory  predicts 
a  larger  bandwidth  for  crosswind  and  this  is  corroborated,  for  instance, 
in  a  comparison  between  Figures  5.11  and  5.12.  The  qualitative  correspondence 
between  theory  and  experiment  is  seen  to  be  perfect  in  these  examples. 
Quantitatively  the  correspondence  between  theory  and  experiment  is 
generally  good  for  not  too  shallow  grazing  angles.  The  less  satisfactory 
bandwidth  comparison  in  Figure  5.13  may  be  caused  by  the  effect  of  shadowing 
as  we  have  already  indicated  in  Chapter  IV,  Figures  4.14. 

In  summary,  we  have  shown  that  the  Kirchhoff  theory  appears  to 
give  quite  accurate  predictions  of  the  frequency  spreading  function  for 
all  surface  roughness  conditions.  We  have  found  that  experimental 
verification  of  the  theory  requires  very  accurate  alignment  of  the 
transducers,  because  even  small  misalignments  which  might  be  caused  by 
transducer  imperfections  can  result  in  significant  spectral  shifts. 

However,  when  all  of  these  factors  are  taken  into  account  the  correspondence 
between  theory  and  experiment  is  seen  to  be  very  good. 


Chapter  VI  Physical  Arguments  for  the  Creation  of  Asymmetric  Sidebands 


6.0  Introduction 

In  the  previous  two  chapters  we  have  used  the  Kirchhoff  integral 
method  in  a  rather  formal  way  to  compute  the  spectral  properties  of  sur- 
f ace-scattered  signal.  We  now  consider  some  physical  arguments  that 
provide  a  more  heuristic  explanation  for  some  of  the  results  obtained 
previously.  We  consider  in  particular  the  two  mechanisms  that  have  been 
shown  to  be  responsible  for  the  spectral  asymmetry  of  the  scattered 
signal. 

6.1  Effect  of  Source-Receiver  Geometry 

Asymmetry  of  the  scattering  geometry,  i.e.  where  the  source  and 
receiver  are  at  different  depths,  has  been  shown  to  result  in  a 
spectral  asymmetry  if  the  source  and  receiver  are  in  an  up  or  downwind 
configuration.  This  asymmetry  has  a  simple  physical  explanation  which 
is  illustrated  in  Figure  6.1,  with  the  source  located  at  a  shallower 
depth  than  the  receiver.  The  water  waves  are  assumed  to  move  from  source 
to  receiver.  For  this  configuration  we  found  earlier  that  the  upper 
doppler  sidebands  have  a  larger  amplitude  than  the  lower  ones. 

Assume  that  the  acoustic  wave  frequency  is  high  enough  so  that 
the  reflecting  surface  can  be  represented  by  a  tangential  plane  at  that 
reflecting  point.  The  planes  at  point  A  and  B  in  Figure  6.1  represent 
two  of  such  surface  tangential  planes,  where  the  reflections  from  these 
planes  will  reach  the  receiver.  Point  A  represents  an  arbitrary  point 
on  the  -Y  axis  and  point  B  represents  the  point  on  the  +Y  axis.  Since 
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Figure  6.1  Asymmetric  downwind  scatter  geometry 

attenuation  factor. 
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the  surface  wave  is  moving  from  source  to  receiver,  the  surface  plane 
at  point  A  moves  downward  and  compresses  the  incoming  acoustic  wave 
and  therefore  produces  an  up-shifted  doppler  frequency  component  in 
the  received  signal  spectrum.  On  the  other  hand,  the  surface  plane 
at  point  B  moves  upward  and  relaxes  the  acoustic  wave  and  produces  the 
down  doppler  shift.  If  the  points  A  and  B  are  chosen  to  produce  equal 
doppler  shifts,  the  upper  (from  point  A)  and  the  lower  (from  point  B) 
doppler  sideband  amplitudes  can  be  obtained  by  the  intensities  of  these 
ray  paths.  Prom  Huy gen's  principle,  each  illuminated  surface  point  be¬ 
comes  the  source  of  new  wave  front,  the  attenuation  factor  for  each 
reflected  ray  is  (Tq^)  1.  A  plot  of  (r^r^)-1  for  all  the  points  on 
the  y  axis  is  shown  above  the  scatter  geometry  in  Figure  6.1.  For  a 
reflection  at  the  specular  point  0,  the  attenuation  factor  is  (rQQr2Q^ 
The  attenuation  factors  for  points  A  and  B  can  be  found  from  the 
correspoinding  points  labelled  A'  and  B'  on  the  (r  r^)  1  curve.  We 
see  that  the  magnitude  of  (r^r^)  ^  at  point  A’  is  larger  than  that 
at  point  B',  therefore  the  upper  sideband  is  larger  than  the  lower 
sideband.  This  conclusion  is  the  same  as  that  reached  in  the  earlier 
chapters. 


The  same  mechanism  operates  in  the  Kirchhoff  integral  which  we 
used  to  develop  the  mathematical  solutions.  The  original  form  of  the 
Kirchhoff  integral  is  given  in  equation  (2.1.6),  repeated  here  for 
convenience 


H(VC)  =  4tt 
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s  3n  {BSBR 
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}•  dS  (6.1.1) 


Observe  that  the  attenuation  factor  is  explicitly  included.  The 
integral  of  equation  (6.1.1)  can  be  read  as  the  sum  of  all  ray  components 


reflected  from  surface  S.  Each  component  is  a  normal  derivative  of  a 

—i  o 

quantity  which  has  amplitude  (r^r^)  and  phase  angle  —  (r^+r^) . 

If  the  direction  of  the  surface  wave  velocity  is  reversed,  the 
plane  at  point  B  will  then  move  downward  and  compress  the  acoustic 
wave.  The  plane  at  point  A  will  relax  the  acoustic  wave  instead  of 
compressing  it.  In  other  words,  the  role  of  plane  A  and  B  are 
switched  and  so  are  the  upper  and  the  lower  sidebands.  Therefore, 
the  direction  of  the  asymmetry  is  reversed. 

A  crosswind  geometry  is  shown  in  Figure  6.2.  On  the  Y  axis, 
two  surface  tangential  planes,  which  give  the  same  amount  of  frequency 
shift  at  the  receiver,  are  shown  at  points  A  and  B.  The  attenuation 
factor  (r^r^)  ^  is  plotted  on  top  of  the  scatter  geometry.  It  is 
observed  that  A  and  B  are  symmetric  and  the  attenuation  factors  at 
these  two  points  are  identical.  In  fact,  the  plot  of  (r^r.^)  *  is 
symmetric  with  respect  to  the  specular  point  0.  Hence  there  should 
be  no  sideband  asymmetries  in  the  crosswind  configuration  if  only  the 
scatter  geometry  is  considered. 

As  indicated  in  the  previous  chapters,  the  upper  and  lower  side¬ 
bands  are  symmetrical  if  Tqq  «*  r^Q  or  if  the  grazing  angle  is  equal  to 
90°.  This  result  is  obvious  since  (rQri)  1  shows  a  symmetric  distri¬ 
bution  under  these  conditions. 

The  reason  why  spectral  asymmetry  has  usually  not  been  observed 
in  scattering  computations  using  the  Kirchhoff  integral  is  that  the 
specular-point  expansion  used  by  Eckart  [31]  and  by  others  [36],  etc. 
replaces  the  factor  (r^r^)  ^  by  a  constant  (r^r^)  ^  for  the  sake  of 
mathematical  convenience.  This  approximation  clearly  removes  the 
asymmetry  contributed  by  the  geometric  factor. 
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Figure  6.2  Asymmetric  crosswind  scatter  geometry  and  attenuation 

factor. 


6.2  Effect  of  the  Surface  Slope  Distribution 

In  the  discussions  of  Chapters  III  and  IV,  we  found  that  slope 
asymmetry  of  the  surface  waves  results  in  asymmetry  in  the  received 
spectrum.  To  show  only  the  effect  of  slopes,  we  assume  a  symmetrical 
scattering  geometry  by  letting  r^  =  r^^  in  an  up  (down)  wind  configura¬ 
tion  shown  in  Figure  6.3.  The  upper  plot  shows  an  assumed  asymmetric 
distribution  of  surface  slopes  for  a  wind  driven  surface.  In  the  lower 
plot  the  tangential  planes  at  A  and  B  are  indicated  as  in  the  previous 
two  figures.  The  arrow  shows  the  direction  of  the  wind  and  of  the 
surface  wave  motion.  The  slope  distributions  corresponding  to  the  two 
surface  points  A  and  B  are  A'  and  B',  respectively.  It  is  seen  the 
slope  distribution  at  point  A'  is  larger  than  that  at  point  B'.  Hence 
there  are  more  positive  slopes  and  fewer  negative  slopes  on  the  surface 
and  we  would  expect  more  energy  in  the  upper  sideband  (reflected  from 
positive  slopes)  than  the  lower  sideband  (reflected  from  the  negative 
slopes) . 

For  the  same  geometry,  if  the  wind  direction  is  reversed,  the 
slopes  at  point  B  will  compress  the  acoustic  wave  but  it  also  will  have 
a  higher  slopes  distribution.  Therefore,  reversing  wind  directions  only 
exchanges  the  roles  of  points  A  and  B,  and  has  no  effect  on  the  direction 
of  the  spectral  asymmetry. 

The  same  result  can  be  obtained  for  the  crosswind  geometry,  which 
is  shown  in  Figure  6.4.  The  reflection  from  surface  plane  A  constitutes 
the  upper  sideband  and  has  more  energy  in  this  band. 
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Figure  6.4  Symmetric  crosswind  scatter  geometry  and  surface 

slope  distribution. 
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When  the  slope  terms  are  retained,  even  approximately,  we  have  shown  that 
slope  asymmetries  result  in  spectral  asymmetries. 
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Chapter  VII  Summary  and  Conclusions 


In  this  study,  we  have  examined  the  problem  of  frequency  spreading 
in  the  acoustic  surface  scattering  from  both  the  mathematical  and  the 
physical  point  of  views.  The  result  of  this  investigation  is  the 
establishment  of  spectral  characteristics  as  a  function  of  three 
effects;  viz.,  scattering  geometry,  surface  statistics,  and  acoustic 
frequency.  Expressions  are  derived  based  on  an  arbitrary  scattering 
geometry  to  allow  for  general  applications. 

The  analysis  considered  two  different  kinds  of  surface  model; 
a  deterministic  one  and  a  random  one.  The  deterministic  model  is 
treated  in  Chapter  III.  A  two-harmonic  surface  Vave  was  used  to  model 
a  wind-driven  surface.  The  spectral  behavior  of  the  acoustic  signal 
scattered  from  this  deterministic  surface  was  derived  for  small  surface 

waveheight.  It  was  found  that  the  surface  scattering  process  introduced 

/ 

doppler  shifted  spectral  lines  into  the  received  signal  spectrum.  The 

I 

amount  of  frequency  shift  is  a  function  of  surface  wave  frequencies. 

The  amplitudes  of  the  two  spectral  lines  adjacent  to  the  transmitted 
acoustic  frequency  are  generally  unequal.  Under  crosswind  condition, 
the  upper  spectral  amplitude  is  always  larger  than  the  lower  one.  In 
the  down  (up)  wind  case,  the  relative  amplitudes  depend  on  wind  direction 
and  relative  depths  of  source  and  receiver.  Scattering  from  a  random 

t 

surface  was  considered  in  Chapter  IV.  The  surface  height  distribution 

was  assumed  to  be  approximately  Gaussian, ibut  such  as  to  permit  a  non- 

1 

zero  third  order  surface  slope  moment.  It  was  found  that  this  moment  is 
one  of  the  two  mechanisms  that  control  the  spectral  amplitudes.  The  other 


mechanism  is  the  surface  scattering-geometry  asymmetry.  A  detailed 
discussion  of  how  the  two  mechanisms  control  the  spectral  amplitudes 
was  given  for  both  the  slightly  rough  and  the  very  rough  surface 
situations.  In  each  situation  we  considered  the  downwind  and  cross- 
wind  cases  separately.  The  result  obtained  from  the  slightly  rough 
random  surfaces  was  compared  against  that  obtained  from  the  deterministic 
model.  We  showed  that  the  latter  could  be  a  limit  case  of  the  former. 

All  these  discussions  are  for  forward  specular  scattering.  In  Chapter 
V,  the  theoretical  and  the  experimental  results  are  examined,  and  the 
amplitude  asymmetries  due  to  misaligned  source-receiver  geometry  (off- 
specular)  are  discussed.  The  major  mathematical  results  of  Chapter 
III,  IV  and  V  were  reviewed  from  a  physical  point  of  view  in  Chapter  VI. 

A  simple  surface  reflection  model  was  used  to  explain  the  analytical 
results  qualitatively. 

It  was  found  that  the  deterministic  model,  random  surface  and 
physical  modelling  all  lead  to  consistent  conclusions  about  the  frequency 
spreading  function.  The  major  results  can  be  summarized  as  follows: 

(a)  The  acoustic  signal  scattered  from  a  sinusoidal  surface  produces 
distinct  spectral  lines  as  sidebands  in  the  received  signal  spectrum. 

The  spectral  lines  are  doppler  shifted  from  the  source  frequency  with 
frequency  shifts  equal  to  the  surface  wave  frequency  and  its  harmonics. 

The  number  of  significant  spectrum  lines  depend  on  the  surface  waveheight 
(the  roughness) .  More  spectral  lines  will  result  if  the  surface  wave- 
height  increases  its  value.  If  the  surface  wave  has  more  than  one 
frequency,  the  resulting  spectrum  will  have  additional  spectral  lines  with 
doppler  shifts  equal  to  the  mixed  frequencies.  The  total  energy  is 
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conserved,  however,  when  more  spectral  lines  are  formed,  the  amplitude 
of  each  spectral  line  decreases. 

9  (b)  When  the  acoustic  signal  Is  scattered  from  a  random  moving  surface, 

the  resulting  spectrum  usually  shows  frequency  spread  around  the  source 
frequency  and  frequency  spreads  around  the  doppler-shif ted  surface  wave 
#  frequencies  (i.e.  the  surface  wave  still  has  one  or  more  characteristic 

frequencies  despite  the  randomness).  The  width  of  the  spread  is  usually 
controlled  by  the  surface  waveheight  and  by  the  surface  slopes  correlation 
1*  constants.  The  smaller  the  correlation  constants  are,  i.e.  the  more 

the  surface  height  or  slopes  are  uncorrelated,  the  larger  is  the 
frequency  spread  (bandwidth).  For  surfaces  with  very  large  roughness, 

**  the  received  signal  shows  a  spectrum  that  has  a  Gaussian  shape.  This 

result  can  be  extended  from  the  small  surface  roughness  situation  by 
widening  the  doppler  shifted  spectrum  as  surface  roughness  increases, 
and  finally  the  bandwidth  of  all  sidebands  become  wide  enough  so  that 
they  merge  with  the  center  spectrum  to  form  a  single  spectrum  around 
the  source  frequency. 

« 

(c)  The  magnitudes  of  the  sidebands  on  both  sides  of  the  acoustic 
frequency  are  in  general  unequal.  Two  different  mechanisms  are 
responsible  for  this  amplitude  difference.  One  of  these  is  the  surface 

* 

slope  asymmetry  which  always  results  in  the  upper  sideband  being  larger 
than  the  lower  sideband.  This  effect  is  independent  of  wind  direction 

or  of  the  relative  depths  of  source  and  receiver.  The  second  mechanism 

t 

is  the  asymmetry  of  the  scattering  geometry;  i.e.  a  difference  in  the 
depths  of  source  and  receiver.  If  the  ratio  between  depths  of 

source  and  receiver  exceeds  2  and  if  the  source-receiver  configura- 

t 

tion  is  down  (or  up)  wind,  the  second  mechanism  becomes  dominant  in 
determining  the  relative  sideband  amplitudes.  Specifically,  the  conditions 

are,  if  the  source  is  at  a  shallower  depth  than  the  receiver  is,  and  if 

I 
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the  wind  blows  from  source  to  receiver  (downwind) ,  then  the  upper 
sideband  is  larger  than  the  lower  sideband.  By  either  reversing 
the  wind  direction  or  maintaining  the  wind  direction  but  placing 
the  source  at  a  larger  depth  than  the  receiver,  the  direction  of 
sideband  asymmetry  is  reversed.  Since  the  second  mechanism  is  a 
wind-direction  dependent  effect,  it  has  no  effect  on  the  sideband 
amplitude  asymmetries  at  crosswind. 

Good  agreement  was  obtained  for  random  surfaces  under  a  variety 
of  roughness  conditions.  Some  of  the  differences  that  were  initially 
found  between  experimental  results  and  the  theory  turned  out  to  be 
the  result  of  small  measurement  inaccuracies.  By  modifying  the 
theory  to  permit  consideration  of  misalignment  in  hydrophone  placement, 
it  was  found  that  very  small  misalignments  could  result  in  fairly 
substantial  doppler  shifts.  Wien  these  misalignments  were  removed, 
excellent  agreement  between  experiment  and  theory  was  obtained  in 
practically  all  cases. 

The  success  of  the  theory  in  matching  experimental  results  can 
be  regarded  as  a  verification  indication  of  the  Kirchhoff  integral 
method  for  dealing  with  the  forward  scattering  problem.  Although  the 
sideband  asymmetries  that  were  one  of  the  major  items  of  interest  in 
this  research  were  not  discovered  by  the  Kirchhoff  method  in  its 
commonly  used  simplified  form,  it  has  been  shown  that  if  some  often 
neglected  factors  were  retained,  the  Kirchhoff  method  gives  experi¬ 
mentally  verifiable  results.  Furthermore,  the  effect  of  surface- 
slope  asymmetries  and  of  small  hydrophone  misalignments  had  not  been 
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previously  predicted.  Thus  although  this  method  has  frequently  been 
criticized,  and  although  it  is  difficult  to  justify  all  the  simplifying 
assumptions  on  which  it  rests,  it  appears  to  be  able  to  give  new  and 
useful  results. 
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APPENDIX  A:  Derivation  of  Helmholtz  Integral 


Consider  a  monochromatic  scalar  wave 


p(x,y,z,t)  =  p  (x,y,z)e“^wt 


(A-l) 


It  satisfies  the  wave  equation 


(y2  +  ^2  )p8<x»y.z>  *  0 


(A- 2) 


Let  p(x,y,z)  be  a  point  within  a  volume  V  bounded  by  a  surface  S 
as  shown  in  Figure  2.1.  If  we  assume  a  function  U(x,y,z)  which  together 
with  its  first  order  derivatives  is  continuous  and  finite  within  and 
on  S,  then  from  Green's  theorem 
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when  n  is  the  surface  normal  defined  positive  outward.  We  also  assume 
U  to  satisfy  the  wave  equation. 


(V2  +  —  )U  =  0  (A-4) 

c 

From  equation  (A-2)  and  (A-4) ,  we  see  that  the  left  side  of 
equation  (A-3)  is  zero  for  every  point  of  V,  therefore  we  have 
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We  now  assume  U  to  be  a  function  of  the  form 

u 
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where  is  the  distance  from  p(x,y,z)  to  any  other  point  in  V. 

We  have  assumed  U  to  be  a  continuous  function  within  V,  and  since 

CO  pk 

i  *  * »  j* 

eJ  c  1  /r^  is  singular  for  r^  *  0,  we  must  exclude  the  point 
p(x,y,z)  from  the  domain  of  integration.  We  surround  P  by  a 
small  sphere  of  radius  e  and  let  the  integration  be  taken  through¬ 
out  the  volume  between  S  and  S',  where  S'  is  the  surface  of  the 
small  sphere  enclosing  P.  We  can  rewrite  equation  (A-5)  as 
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Hence 
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Notice  that  the  first  term  and  the  third  term  in  the  kernel  of  the 

right-hand  side  integral  are  proportional  to  ,  where  S'  is 
2 

proportional  to  e  ,  therefore  if  we  let  s  0,  the  first  and 
the  third  term  in  the  right-hand  side  integral  gives  no  contribution 
to  the  integrated  value.  However,  the  integration  of  the  second 

1  X 

term  yields  a  value  4xpg,  since  eJ  c  closes  to  unity  when  e  0. 
Dividing  both  sides  of  equation  (A-8)  by  4ir  gives  equation  (2.1.1) 
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APPENDIX  B:  Surface  Integral  Transformation 


Consider  a  surface  integral  I  of  the  form 


I 


f (x,y,S(x,y))* 


ds 


(B-l) 


where  S  represents  the  surface  area  over  which  the  integral  is 
to  be  performed,  d£  is  a  small  surface  element  whose  direction 
is  that  of  the  surface  normal;  x,  y  and  S (x,y)  are  the  coordinates 
of  an  arbitrary  point  on  the  surface  S  in  an  x,y,z  coordinate 
system  as  shown  in  Figure  B.l.  The  location  of  each  point  of 
the  surface  area  S  is  completely  specified  by  the  vector 


£(x,y)  =  xx+y^  +  C  (x,y)z 


(B-2) 


Let  S  be  a  smooth  surface  and  assume  r(x,y)  to  have 
continuous  first  derivatives  along  the  x  and  y  directions. 
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These  two  partial  derivatives  define  the  slope  of  the  tangential 
plane  of  S  at  the  surface  point  (x,y (x,y)) .  Therefore  the  surface 
normal  n  at  point  (x,y,C(x,y))  is 
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The  angle  0  between  surface  normal  n  and  the  z  axis  is  the  arc 
cosine  of  the  inner  product  of  n  and  Is,  or 
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If  the  surface  element  is  sufficiently  small,  we  can  replace 
|  ds|  by  dxdy/cos  9;  hence 
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The  Integral  in  equation  (B-l)  can  now  be  rewritten  as 

ffs  [k*  +  h*  +  h*  £^T0~  (B"8) 
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where  S  represents  the  projection  of  s  on  the  x-y  plane  as  shown 

xy 

g 

in  Figure  B.l.  - —  has  been  replaced  by  its  equivalent  expression 
on 

in  x,y,z  coordinates. 
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APPENDIX  C:  Spectrum  Computation-Deterministic  Surface,  Downin 


From  equation  (3.2.7)  and  C  =0,  H(m , t)  can  be 

x  Q 


rewritten  as 
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e  c  dxdy 
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A  =  — — 


_ j  “0 

L2(d'0)  cr 
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F  =  2  siniji 


(C-2) 


Since  C(x,y,t)  is  not  a  function  of  x,  the  integral  with 

* 

respect  to  x  in  equation  (C-l)  can  be  easily  performed  and  H(a)Q,t) 
becomes 


ja)oTs  r  u 
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(C-3) 


By  use  of  the  one  dimensional  surface  wave  model  in  equation 
(3.1.1),  ^(x,y,t)  can  be  put  into  the  equivalent  form 
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4(x,y,t)  -  -±  (ejpy-jQt+J*  +  e'Jpy+jfit-^) 
2 


+  -§  (eJ2py-J2ftt  +  e“j2py+j2nt  }  (c_4) 


The  slope  along  the  y  direction  is  obtained  by  differentiation 
of  equation  (C-4) ,  and  it  becomes 
ph. 


C/x.y.t,  -  i  V 


+  jph2(e)2py-32nt  -  (C-5) 


and  therefore 
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(C-6) 


The  waveheight  of  the  water  surface  is  assumed  to  be 
small  compared  to  the  acoustic  wavelength.  From  equation  (3.2.8) 


e  Jc  F  ~  !  +  j  "o  F  b:[ejpy-jfit+j^+  e-jpy+jfit-j((), 
c  2 1 
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Substituting  equations  (C-4)  to  (C-7)  into  equation  (C-3) 
results  in 
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Since  HC0EF  has  the  form  cQ  +  c^,  the  integral  in 
equation  (C-8)  can  be  integrated  in  closed  form.  The  result  of 
this  integration  is 
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APPENDIX  D:  Proof  of  X  >  0 

The  water  surface  function  is 

4(x,y,t)  ■  cos(py-ftt+$)+h2  cos(2py-2flt)  (D-l) 

which  has  a  slope  of  the  form 


y  3y 


=  -  ph^  sin(py-Qt+$)  -  2ph^  sin(2py-2Rt)  (D-2) 


since  the  water  wave  is  assumed  to  be  one-dimensional  in  the  y 
direction,  the  slope  is  zero  in  x  direction. 


4  is  given  by 


4y2  =  p2^2  sin2(py-flt+<j>)+4p2h22sin2(2py-2nt) 


+  4P  h^h28in(py-£2t+^)sin(2py-2fit)  (D-3) 


2tt 

In  performing  the  average  over  T  «  ,  the  cross-product 

term  vanishes .  Hence 


2U  2 

2  P  "l  2  2 

V  "  ~~T~  +  2p  V 


(D-4) 


For  hydrodynamic  reasons,  the  water  surface  slope  is  in 
general  less  than  unity  in  absolute  magnitude;  i.e.  slope  angles 
are  generally  less  than  45°.  Therefore 
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0  <  — ~  +  2pZh2  <  1 
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use  of  this  inequality  in  equation  (3.2.14),  viz. 
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results  In 


>  hi  .  vln*o  r  p2hi2 

+  r„„  c  L  4 


00  ‘10 


Since  the  grazing  angle  iJIq  is  between  0  and 
therefore  X  >  0. 
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APPENDIX  E:  Spectrum  computation  -  Deterministic  Surface,  Crosswind 


The  distances  r^  and  r^  can  be  expressed  in  the  following 


ro  =  J  (x"r00CO8V2  +  y2  +  (t+  roo  ain  V2 


(E-l) 


ri  “  J +  y2  +  (?+r10sinij>0)2 


(E-2) 


The  Fresnel  expansion  of  equation  (E-l)  and  (E-2)  gives,  cf. 
equations  (2.2.6) 


ro  “  roo  ■  x  cos  ^0  +c  sin  *0  + 


2.2  .  2 

x  sin  +  y 
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r^  *  r^Q  +  x  cos  iJ/q  +C  sin  t|»Q  + 


2  2  2 
x  sin  V'q  +  y 


(E-4) 


and  the  total  transmission  distance  r^  +  r^  is  the  sum  of  these 


two  equations 


ro  +  ri  “  roo  +  rio  +  25  8i“*o  + 


2.2  .  2 

x  sin  t(»Q  +  y 


(E-5) 


2r  r 

where  R  ■  - — -  was  defined  in  equation  (2.2.14). 

r00+r10 

The  attenuation  factor  — —  can  be  approximated  by 

0  1 

(see  equation  (2.2.11)) 

77~  -  exp |  (  —  -  ■—)  x  cos  *Q|  (E-6) 

r0rl  r00  10  '  00  10  ’ 

According  to  equation  (3.2.6),  the  beam  pattern  function  on  the 

x-y  plane  can  be  written  as 
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bsbr  -  exP 


{  ■  ~2~^ -  l 

1  L  K) 


x2sin2^0  +  y2 ]| 


(E-7) 


Substituting  equations  (E-5)  to  (E-7)  into  the  Kirchhoff  integral 
formula,  one  gets 

-Bx-Ay-Ex+j^%; 


H(uQ,t) 


juJ-T  /?  r<*> 

trr~  J  l1?F(1-!>'rsli 

00  10  J-aJ-a> 


dxdy 


(E-8) 


where  A,  B,  E  and  F  are  defined  in  Appendix  C. 


Comparisons  of  equations  (C-l)  and  (E-8)  shows  that  these  two 
integrals  are  very  similar  in  form  except  that  in  equation  (E-8) , 
there  is  an  e  term  in  place  of  e  y  ,  also  2By  +  E  is  replaced  by 
2Ay.  Integrating  equation  (E-8)  with  respect  to  x  gives 
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4^rB 


roorio 


F(l-ty2)+2 


2  0 

-Ay/+j  ~  f; 


dy  (E-9) 


Using  the  surface  wave  model  of  equation  (3.1.1),  we  can 
evaluate  the  integral  in  equation  (E-9)  in  the  same  way  as  in  Appendix 
C.  The  result  is 
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H(o>0»t) 


where 


4  ^B 


-  {Ho  +  \  [ve"J"t+ 

roorio  '  "-1 


(E-10) 


2.  2 
P  h. 


2  ,  2.  2. 
ii>n  „  3p  h.  h. 


*0-  J  ^  F(1  -  .  2p2h22)+  f2  _T^2  cos2^ 

2  c  4 


E-2 


r 


APPENDIX  Fj  Representation  of  Transfer  Function  H(u)g,t)  in  a 

General  Scatter  Geometry 

The  form  of  the  channel  transfer  function  is 


B(“o,t)  ■  -sUL  {  k  -  h  -  S  U  {BA  *'  ~-ri-~>dy  (F-1) 


where  Tq  and  can  be  expressed  as 
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(F-2) 


-x  “  J(.x-t1qcob  4iRsin$R) 2+  (y+r x ^co Rcos  $  R) "^+(^+i~s (F-?) 


As  in  Chapter  II,  the  Fresnel  expansion  of  r^  and  r^  have  the 
forms 
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Therefore  Tq+t^  can  be  obtained  from  the  sum  of  these  two  equations 
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The  attenuation  factor  -  is  approximated  as  (see  Chapter  II) 

rori 

l  1  cosi(>Tsin<f>T  cos\J>Rsin<t>R  cos^cos^  cosi|>rcos4>r 
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(P-7) 


The  beam  pattern  function  on  the  x-y  plane  can  be  computed  as  in 
Chapter  II,  equation  (2.3.13) 


Mr  "  exp{ - 2 - ? - f22 —  i 

'■J <V  vlV 


(F-8) 


where  L  (<0 ,  L  (wn)  and  L  (un)are  defined  as  in  equation  (2.3.14). 
xx  0  yy  0  xyu  0 
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If  we  express  B„B  -  in  the  general  form 
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APPENDIX  G:  Sign  of  Third-Order  Slope  Moment 


In  general,  the  surface  slope  on  the  windward  side  is 
shallower  than  that  on  the  leeward  side.  To  demonstrate  that 
the  third-order  slope  moment  is  always  negative,  we  consider 
an  idealized  surface  profile  as  shown  in  Figure  G.l. 


The  wind  direction  is  from  left  to  right.  The  slope  on 
the  windward  side  is  |  BD/AD  j  and  on  the  leeward  side  it  is 
-  j  BD/CD  j  .  The  underline  represents  the  distance  between 
the  two  points.  By  hypothesis 
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so  that 

|  AD  |>|CD|  (G-2) 

The  third  order  moment  is  defined  as  the  average  of  the 
cube  of  the  surface  slope.  The  idealized  surface  wave  is  periodic. 


therefore  can  be  calculated  as 
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APPENDIX  H:  Power  Spectrum-Random  Surface,  Slightly  Rough 


We  first  rewrite  J(z  ,r  )  in  equation  (A. 1.7)  by 

yl  y2 

computing  the  average  of  the  product  of  two  brackets. 
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Make  a  change  of  variable  x2  =  x^  +  £,  the  integration  over  x^  and 
C  in  equation  (A. 1.6)  is  easily  performed. 
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We  make  another  change  of  variable 
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and  define 
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Performing  the  integration  over  y.^  results  in 
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If  we  now  substitute  equations  (4.1.17)  and  (4.2.1)  into  equation 
(H-7),  $(o)0>t)  can  be  written  as 
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The  remaining  integration  over  n  are  straightforward,  but 
involve  a  lot  of  algebra.  The  results  are 
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The  final  signal  spectrum  is 
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Again,  we  have  defined  new  symbols  to  simplify  the  expression, 
they  are 
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Hence  r(u,co^) 


can  finally  be  expressed  as  a  three-part  spectrum 


r(u>,ioo)  =  K1[Spo(u,o)o)+  Spi+(u,iO+  Sp^_(u),u)^)  ] 


(H-14) 


I 


# 


H-7 


t2V*0)2 


SPO(“,“o)  =  f006(h>-W0)+^*T2e2f20e 


T3£3f30e 


t32(^0)2 


(H-15a) 


SPl+(“ 


•"o’  ■{' 


Tq2  (w-oig+PCp) 2 


J-8-  /if  t  f  e 
<  2  0  0+ 


^  T0TO 

— i‘2V 


(  T02+T22)! 


l02T22(ia-ypCp)2 

4(T02  +  T22) 

I02T3Z<u-»(jfpCp)2 


.  *  T0T3 

+  - - - -  £,f-.e 

Vi  3  3± 


4(T02+T32) 


(T02+T32  )! 


(H-15b) 


Case  I.  Down-Wind  Geometry 
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and  the  grazing  angles  of  transmitter  and  receiver  are  the  same 
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If  the  beamwidth  L(o^)  is  larger  than  the  dimension  of  the 
Fresnel  zone,  then  from  [75] 
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Note  that  M  ■  E  is  real,  therefore  from  equation  (H-13) 
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Case  II  Cross-Wind 
*T  =  90° 

^R  =  270° 
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All  the  assumptions  and  variables  are  identical  to  those 
of  Case  I,  except  we  now  have 
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APPENDIX  I:  Power  Spectrum-Random  Rough  Surfaces 


The  integration  over  x^,  x2  and  are  the  same  as  in 
Appendix  H.  The  result  is  equation  (H-7) 
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Carrying  out  the  integration  over  n  and  applying  equation  (H-ll) , 
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where  d  has  the  same  definition  as  in  equation  (H-13)  and 
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The  signal  spectrum  is 
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Case  I .  Downwind 

The  conditions  in  equations  (H-16~20)  are  still  valid.  Let's 
consider  a  large  beam  situation  where 
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From  equation  (1-3) 
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Therefore  from  equation  (1-4) 


Case  II.  Crosswind 


Combining  equations  (H-25-  27)  and  the  conditions  in  the  down¬ 
wind  geometry,  the  signal  spectrum  of  cross-wind  geometry  is 
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APPENDIX  J:  Signal  Processing  Technique 


The  acoustic  source  signal  is  shown  on  the  left  side  of 
Figure  J.l.  It  consists  of  a  series  of  narrow  pulses,  with  a  time 
interval  of  10  msec  between  pulses.  The  transmitted  signal  is 
wideband  with  bandwidth  ranging  from  0.2  MHz  to  1.2  MHz. 

The  signal  received  at  the  receiver  is  also  in  pulse  form 
corresponding  to  each  pulse  sent  at  the  source.  The  received 
signal  h(t,t)  is  sampled  and  recorded.  Each  received  pulse  is 
sampled  at  a  sampling  interval  of  350  nsec,  256  samples  per  pulse. 

The  beginning  time  of  each  received  pulse  is  recorded  in  t^,  with  m 
ranging  from  1  to  64.  A  data  file  consists  of  64  x  256  =  16,384 
samples . 

The  signal  spectrum  computational  technique  is  shown  in  Figure 
J.2.  Each  data  file  is  first  passed  through  a  256  point  FFT  for 
each  t^.  For  350  ns  sampling  time,  this  gives  a  frequency  range 
from  0  Hz  to  1.42  MHz  with  a  frequency  resolution  of  11.16  KHz. 

One  of  the  256  FFT  outputs  corresponding  to  the  desired  acoustic 
frequency  u>q  is  selected  for  each  tm.  The  64  values  of  H(a)g,t)  in 
each  file  are  passed  through  a  64  point  FFT,  the  output  of  which  is 
the  raw  signal  spectrum  for  that  file.  This  has  a  frequency  resolution 
of  1.5  Hz.  The  final  spectrum  r(w,u>g)  is  an  average  of  the  outputs  of 
N  of  such  files. 


J-l 


WINO  DIRECTION 


Figure  J.l  Signal  representation  of  acoustic  signal  transmitted  and  received  underwater. 


Figure  J.2  Signal  spectrum  computing  block  diagram 


REFERENCES 


[1]  F.M.  Labianca,  and  E.Y.  Harper,  "Sideband  Structure  of  Sound 
from  a  Harmonic  Point  Source  Scattered  by  a  Rough  Surface 
Moving  over  an  Upward-refracting  Ocean",  J.  Acoust.  Soc.  Am. 

Vol.  61,  No.  2,  378-389,  1977. 

[2]  W.A.  Kuperman,  "On  Scattering  from  a  Moving  Rough  Surface 
Without  a  Kirchhoff  or  a  Farfield  Approximation",  J.  Acoust. 

Soc.  Am.,  Vol.  64,  No.  4,  1113-1116,  1978. 

[3]  J.G.  Zornig,  "Frequency  Spreading  in  Forward  Surface  Scattering", 
Tech.  Report  #7710,  Dept.  Engineering  and  Applied  Science,  Yale 
University,  New  Haven,  Ct.,  December  1977. 

[4]  B.E.  Parkins,  "Scattering  from  the  Time-Varying  Surface  of  the 
Ocean",  J.  Acoust.  Soc.  Am.,  Vol.  42,  1262-1267,  1967. 

[5]  L.  Scharf  and  R.  Swarts,  "Acoustic  Scattering  from  a  Stochastic 
Sea  Surface",  J.  Acoust.  Soc.  Am.,  Vol.  55,  247-253,  1974. 

[6]  C.W.  Horton  and  T.G.  Muir,  "Theoretical  Studies  on  the  Scattering 
of  Acoustical  Waves  from  a  Rough  Surface",  J.  Acoust.  Soc.  Am., 
Vol.  41,  627-634,  1967. 

[7]  C.S.  Clay  and  H.  Medwin,  "Dependence  of  Spatial  and  Temporal 
Correlation  of  Forward-Scattered  Underwater  Sound  on  the 
Surface  Statistics.  I.  Theory”,  J.  Acoust.  Soc.  Am.  Vol.  47, 
1412-1418,  1970. 

[8]  J.F.  McDonald  and  F.B.  Tuteur,  "Moment  Characterization  of  a 
Doubly  Spread  Surface  Scatter  Channel  at  High  Rayleigh  Parameters" 
IEEE  Proc.,  Vol.  62,  No.  11,  1606-1608,  Nov.  1974. 

[9]  l.M.  Fuks,  "Spectral  Width  of  Signals  Scattered  by  a  Disturbed  Sea 
Surface",  Sov.  Phys.  Acoust.,  Vol.  20,  276-281,  1974. 

[10]  J.F.  McDonald  and  P.M.  Schulthelss,  "Asymptotic  Frequency  Spread 
in  Surface-Scatter  Channels  at  Large  Rayleigh  Numbers",  J.  Acoust. 
Soc.  Am.,  Vol.  57,  No.  1,  160-164,  1975. 

[11]  R.L.  Swarts  and  C.J.  Eggen,  "Simplified  Model  of  the  Spectral 
Characteristics  of  High-Frequency  Surface  Scatter",  J.  Acoust. 

Soc.  Am.,  Vol.  59,  No.  4,  846-851,  April  1976. 

[12]  F.B.  Tuteur,  J.G.  Zornig  and  H.  Tung,  "Scattering  of  Underwater 
Acoustic  Signals  from  a  Rough,  Moving  Surface",  Final  Report 
under  Contract  No.  N00014-75-C-0298,  S&IS  Report  No.  7705,  June 
1977. 

[13]  C.  Gazanhes,  J.  Leandre  and  J.P.  Lefebvre,  "Spectral  Structure  of 
an  Ultrasound  Wave  Scattered  by  a  Random  Surface:  Application  to 
the  Scattering  of  Sound  from  the  Sea  Surface",  J.  Acoust.  Soc.  Am. 
Vol.  63,  No.  5,  1347-1352,  1978. 


I 


R-l 


[14]  C.S.  Clay  and  H.  Medwln,  "High-Frequency  Acoustical  Reverberation 
from  a  Rough-Sea  Surface",  J.  Acoust.  Soc.  Am.,  Vol.  36,  2131- 
2134,  1964. 

[15]  D.  Middleton,  "Doppler  Effects  for  Randomly  Moving  Scatterers 
and  Platforms",  J.  Acoust.  Soc.  Am.,  Vol.  61,  1231-1250,  1977. 

[16]  W.I.  Roderick  and  B.F.  Cron,  "Frequency  Spectra  of  Forward- 
Scattered  Sound  from  the  Ocean  Surface",  J. Acoust.  Soc.  Am.,  Vol. 

48,  No.  3,  Pt.  2,  759-766,  Sept.  1970. 

[17]  M.V.  Brown  and  G.V.  Frisk,  "Frequency  Smearing  of  Sound  Forward- 
Scattered  from  the  Ocean  Surface",  J.  Acoust.  Soc.  Am.,  Vol.  55, 
744-749,  1974. 

[18]  J, A.  Shooter  and  S.K.  Mitchell,  "Observations  of  Acoustic  Side¬ 
bands  in  CW  Tones  Received  at  Long  Ranges",  J.  Acoust.  Soc.  Am., 

Vol.  60,  829-832,  1976. 

[19]  l.P.  Lysanov,  "Theory  of  the  Scattering  of  Waves  at  Periodically 
Uneven  Surfaces",  Sov.  Phys.  Acoust.,  Vol.  4,  No.  1,  1-10,  1958. 

[20]  L.  Fortuin,  "Survey  of  Literature  on  Reflection  and  Scattering 
of  Sound  Waves  at  the  Sea  Surface",  J.  Acoust.  Soc.  Am.,  Vol. 

47,  No.  5,  Pt.  2,  1970. 

[21]  Lord  Rayleigh,  The  Theory  of  Sound,  (Dover,  New  York,  1945). 

[22]  J.L.  Uretsky,  "Reflection  of  a  Plane  Sound  Wave  from  a  Sinusoidal 
Surface",  J.  Acoust.  Soc.  Am.,  Vol.  35,  1293-1294  (L) ,  1963. 

[23]  H.W.  Marsh,  "In  Defense  of  Rayleigh's  Scattering  from  Corrugated 
Surfaces",  J.  Acoust.  Soc.  Am.,  Vol.  35,  1835  (L) ,  1963. 

[24]  R.F.  Millar,  "On  the  Rayleigh  Assumption  in  Scattering  by  a 
Periodic  Surface",  Proc.  Cambridge  Phil.  Soc.,  Vol.  65,  773- 
791,  1969. 

[25]  J.G.  Parker,  "Reflection  of  Plane  Sound  Waves  from  an  Irregular 
Surface",  J.  Acoust.  Soc.  Am.,  Vol.  28,  672-680,  1956. 

[26]  S.R.  Murphy  and  G.E.  Lord,  "Scattering  from  a  Sinusoidal  Surface — 

A  Direct  Comparison  of  the  Results  of  Marsh  and  Uretsky",  J.  Acoust. 
Soc.  Am.,  Vol.  36,  1598-1599  (L),  1964. 

[27]  G.R.  Barnard  et  al.,  "Underwater-Sound  Reflection  from  a  Pressure- 
Release  Sinusoidal  Surface",  J.  Acoust.  Soc.  Am.,  Vol.  39,  1162-1169, 
1966. 

[28]  W.C.  Meecham,  "Variational  Method  for  the  Calculation  of  the  Dis¬ 
tribution  of  Energy  Reflected  from  a  Periodic  Surface.  I.",  J. 

Appl.  Phys.,  Vol.  27,  No.  4,  361-367,  1956. 


R-2 


[29]  W.C.  Meecham,  "Fourier  Transform  Method  for  the  Treatment  of 
the  Problem  of  the  Reflection  of  Radiation  from  Irregular 
Surfaces",  J.  Acoust.  Soc.  Am.,  Vol.  28,  No.  4,  370-377,  1956. 

[30]  H.S.  Heaps,  "Reflection  of  Plane  Waves  of  Sound  from  a  Sinusoidal 
Surface”,  J.  Appl.  Phys.,  Vol.  28,  815-818,  1957. 

[31]  C.  Eckart,  "The  Scattering  of  Sound  from  the  Sea  Surface",  J. 
Acoust.  Soc.  Am.,  Vol.  25,  566-570,  1953. 

[32]  B.B.  Baker  and  E.T.  Copson,  The  Mathematical  Theory  of  Huygen's 
Principle,  (Oxford  University  Press,  New  York  1950),  Second 
Edition,  Chapter  2. 

[33]  W.C.  Meecham  "On  the  Use  of  the  Kirchhoff  Approximation  for  the 
Solution  of  Reflection  Problems",  J.  Rational  Mech.  Anal.,  Vol. 

5,  323-333,  1956. 

[34]  D.  Mintzer,  "Discussion  of  the  Paper  by  C.  Eckart  on  Sea  Surface 
Scattering",  J.  Acoust.  Soc.  Am.,  Vol.  25,  1015  (L),  1953. 

[35]  J.  Feinstein,  "Some  Stochastic  Problems  in  Wave  Propagation.  I", 
Trans.  IRE  AP-2,  23-30,  1954. 

[36]  E.P.  Gulin,  "Amplitude  and  Phase  Fluctuations  of  a  Sound  Wave 
Reflected  from  a  Sinusoidal  Surface",  Sov.  Phys.— Acoust. ,  Vol. 

8,  223-227,  1963. 

[37]  D.R.  Melton  and  Sr.  C.W.  Horton,  "Importance  of  the  Fresnel 
Corrections  in  Scattering  from  a  Rough  Surface'.'  I.  Phase  and 
Amplitude  Fluctuations",  J.  Acoust.  Soc.  Am.,  Vol.  47,  290- 
298,  1970. 

[38]  J.F.  McDonald,  "Fresnel-Corrected  Second-Order  Inter frequency 
Correlations  for  a  Surface-Scatter  Channel",  IEEE  Trans.  Commu., 
Vol.  COM-22,  No.  2,  138-145,  1974. 

[39]  R.L.  Swarts  and  L.L.  Scharf,  "Note  on  Acoustic  Scattering  from 
Random-Amplitude/Random-Phase  Sinusoidal  Sea  Surfaces",  J. 

Acoust.  Soc.  Am.,  Vol. 56,  No.  2,  686-687,  1974. 

[40]  F.B.  Tuteur,  J.F.  McDonald  and  H.  Tung,  "Second-Order  Statistical 
Moments  of  a  Surface  Scatter  Channel  with  Multiple  Wave  Direction 
and  Dispersion,"  IEEE  Trans.  Commu.,  Vol.  COM-24,  No.  8,  820-831, 
1976. 

[41]  M. A.  Isakovitch,  "Scattering  of  Acoustic  Waves  at  Small  Inhomo¬ 
geneities  in  a  Waveduct",  Sov.  Phys . -Acoust . ,  Vol.  3,  35,  1957. 

[42]  V.G.  Bezrodnyi  and  I.M.  Fuks,  "Amplitude  and  Phase  Fluctuations 
of  the  Field  in  a  Waveguide  Having  a  Rough  Wall",  Sov.  Phys.- 
Acoust.,  Vol.  17,  No.  4,  444-449,  1972. 


R-3 


[43]  J.R.  Wait,  "Perturbation  Analysis  for  Reflection  from  Two- 
Dimensional  Periodical  Sea  Waves",  Radio  Sci.,  Vol.  6,  387- 
391,  1971. 

[44]  E.Y.  Harper  and  F.M.  Labianca,  "Perturbation  Theory  for 
Scattering  from  a  Point  Source  by  a  Moving  Rough  Surface 
in  the  Presence  of  Refraction",  J.  Acoust.  Soc.  Am.  Vol. 

57,  No.  5,  1044-1051,  1975. 

[45]  F.M.  Labianca  and  E.Y.  Harper,  "Connection  Between  Various 
Small-waveheight  Solutions  of  the  Problem  of  Scattering  from 
from  the  Ocean  Surface",  J.  Acoust.  Soc.  Am.,  Vol.  62,  1144- 
1157,  1977. 

[46]  W.A.  Kuperman,  "Coherent  Component  of  Specular  Reflection  and 
Transmission  at  a  Randomly  Rough  Two-fluid  Interface",  J. 

Acoust.  Soc.  Am.,  Vol.  58,  365-370,  1975. 

[47]  D.  Middleton,  "A  Statistical  Theory  of  Reverberation  and 
Similar  First-Order  Scattered  Fields-Part  Is  Waveforms  and 
the  General  Process,"  Trans.  IEEE  Information  Theory,  Vol. 

13,  372-392,  1967. 

[48]  W.J.  Jobst,  "An  Application  of  Poisson  Process  Models 

to  Multipath  Sound  Propagation  of  Sinusoidal  Signals",  J.  Acoust. 
Soc.  Am.,  Vol.  57,  No.  6,  Pt.  2,  1409-1412,  June  1975. 

[49]  P.  Beckmann,  "Scattering  by  Non-Caussian  Surfaces",  Trans.  IEEE 
Antennas  Propaga. ,  Vol.  AP-21,  No.  2,  169-175,  1973. 

[50]  J.F.  McDonald  and  R.C.  Spindel,  "Implications  of  Fresnel 
Correlations  in  a  Non-Gaussian  Surface  Scatter  Channel", 

J.  Acoust.  Soc.  Am.,  Vol.  50,  No.  3,  Part  1,  746-757,  1971. 

[51]  B.  Kinsman,  Wind  Waves,  Their  Generation  and  Propagation  on 

the  Ocean  Surface,  (Prentice-Hall,  Englewood  Cliffs,  N.J.,  1965) 

[52]  D.E.  Weissman,  "Two  Frequency  Radar  Interferometry  Applied  to 
the  Measurement  of  Ocean  Wave  Height",  Trans.  IEEE  Antennas 
Propaga.,  Vol.  AP-21,  649-656,  1973. 

[53]  J.G.  Zornig,  "Measurement  of  the  Statistics  of  a  Surface  Scatter 
Channel”,  Ph.D.  Dissertation,  Dept,  of  Engineering  and  Applied 
Science,  Yale  University,  New  Haven,  Ct.,  1974. 

[54]  R.C.  Spindel,  "An  Experimental  Investigation  of  a  Statistically 
Distributed  Scattering  Surface  and  Acoustic  Scattering",  Ph.D. 
Dissertation,  Dept,  of  Engineering  and  Applied  Science,  Yale 
University,  New  Haven,  Ct.,  1971. 

[55]  C.  Cox,  and  W.  Munk,  "Measurement  of  the  Roughness  of  the  Sea 
Surface  from  Photographs  of  the  Sun's  Glitter",  J.  Opt.Soc. 

Am.,  Vol.  44,  838-850,  1954. 


R-4 


[56]  A.H.  Schooley,  "A  Simple  Optical  Method  for  Measuring  the 
Statistical  Distribution  of  Water  Surface  Slopes",  J.  Opt. 

Soc.  Am.,  Vol.  44,  37-40,  1954. 

[57]  J.  Wu,  "Directional  Slope  and  Curvature  Distributions  of 
Wind  Waves",  J.  Fluid  Mech. ,  Vol.  79,  Pt.  3,  463-480,  1977. 

[58]  J.G.  Zornig  and  J.  Snyder,  unpublished  results. 

[59]  P.  Beckmann,  "Shadowing  of  Random  Rough  Surfaces",  Trans. 

IEEE  Antennas  Propaga,,  AP-13,  384-388,  1965. 

[60]  P.  Beckmann,  "Scattering  by  Composite  Rough  Surfaces", 

IEEE  Proc. ,  1012-1015,  1965. 

[61]  R. A.  Brockelman  and  T.  Hagfors,  "Note  on  The  Effect  of 
Shadowing  on  The  Backscattering  of  Waves  from  Random  Rough 
Surface",  IEEE  Trans.  Antennas  Propag. ,  AP-14,  621-629,  1966. 

[62]  R.J.  Wagner,  "Shadowing  of  Randomly  Rough  Surfaces",  J. 

Acoust.  Soc.,  Am.  Vol.  41,  138-147,  1967. 

[63]  P.J.  Lynch  and  R.J.  Wagner,  "Energy  Conservation  for  Rough- 
Surface  Scattering",  J.  Acoust.  Soc.  Am.,  Vol.  47,  816-821, 

1970. 

[64]  P.J.  Lynch,  "Curvature  Corrections  to  Rough-Surface  Scattering 
at  High  Frequencies",  J.  Acoust.  Soc.  Am.,  Vol.  47,  No.  3, 

Part  2,  804-815,  1970. 

[65]  B.G.  Smith,  "Geometrical  Shadowing  of  a  Random  Rough  Surface", 
Trans.  IEEE  Antennas  Propaga.,  AP-15,  668-671,  1967. 

[66]  R.R.  Gardner,  "Acoustic  Backscattering  from  a  Rough  Surface 
at  Extremely  Low  Grazing  Angles",  J.  Acoust.  Soc.  Am.,  Vol. 

53,  848-857,  1973. 

[67]  J.C.  Novarini,  and  H.  Medwin,  "Diffraction,  Reflection  and 
Interference  during  Near-grazing  and  Near-normal  Ocean  Sur¬ 
face  Backscattering",  J,  Acoust.  Soc.  Am.,  Vol.  64,  No.  1, 
260-268,  July  1978. 

[68]  V.P.  Glotov  and  Y.P.  Lysanov,  "Relative  Contribution  of  Sur¬ 
face  Mr  Bubbles  and  Waves  to  the  Formation  of  Sea  Reverberation" 
Sov.  Phys.  — Acoust.,  Vol.  14,  311-314,  1969. 

[69]  I.  Tolstoy  and  C.S.  Clay,  Ocean  Acoustics,  (McGraw-Hill, 

New  York,  1966) . 

[70]  A.H.  Schooley, "Simple  Tools  for  Measuring  Wind  Fields  above 
Wind-Generated  Water  Waves",  J.  Geophys.  Res.,  Vol.  68,  No. 

19,  5497-5504,  1963. 


[71]  M.  Born  and  E.  Wolf,  Principles  of  Optica  (Pergamon, 

New  York,  1970).  Chap.  8,  pp.  448-452. 

[72]  L.L.  Beranek,  Acoustics,  (McGraw-Hill,  New  York,  1954). 

[73]  R.J.  Urick,  Principles  of  Underwater  Sound  for  Engineers, 
(McGraw-Hill,  New  York,  1967). 

[74]  J.F,  McDonald,  "Passive  Detection  and  Tracking  Using 
Surface  Scattered  Signals",  Ph.D.  Dissertation,  Dept,  of 
Engineering  and  Applied  Science,  Yale  University,  New  Haven, 
Conn.  1970. 

[75]  P.  Beckmann  and  A.  Splzzichino,  The  Scattering  of  Electro¬ 
magnetic  Waves  from  Rough  Surfaces,  (MacMillan,  New  York, 

1963) . 

[76]  J.G.  Parker.  "Reflection  of  Plane  Sound  Waves  from  a 
Sinusoidal  Surface",  J.  Acoust.  Soc.  Am.,  Vol.  29,  377- 
380,  1957. 

[77]  W.C.  Meecham,  "Point  Source  Transmission  Through  a  Sinusoidal 
Ocean  Surface",  J.  Acoust.  Soc.  Am.,  Vol.  64,  No.  5,  1478- 
1481,  1978. 

[78]  E.O.  LaCasce  and  P.  Tamarkin,  "Underwater  Sound  Reflection 
from  a  Corrugated  Surface",  J.  Appl.  Phys.,  Vol.  27,  138- 
148,  1956. 

[79]  J.G.  Zornig  and  J.F.  McDonald,  "Direct  Measurement  of  Surface- 
Scatter  Channel  Coherence  by  Impulse  Probing",  J.  Acoust.  Soc. 
Am.,  Vol.  55,  1205-1211,  1974. 

[80]  H.  Schwarze,  "Model  for  the  Frequency  Spread  of  Backscattered 
Underwater  Sound  Based  on  The  Facet  Model  of  The  Rough  Sea 
Surface,"  J.  Acoust.  Soc.  Am.,  Vol.  64,  No.  2,  605-613,  Aug. 
1978. 

[81]  P.M.  Sc hultheiss,  J.  Snyder,  H.  Tung  and  J.G.  Zornig,  "Acoustic 
Scattering  from  a  Wind  Driven  Surface",  Final  Report  under 
Contract  N00014-75-C-1014,  Dept,  of  Engineering  and  Applied 
Science,  Yale  University,  June  1976. 

[82]  W.B.  Davenport  and  W.L.  Root,  An  Introduction  to  the  Theory 
of  Random  Signals  and  Noise,  (McGraw-Hill,  New  York,  1958) . 

[83]  S.0.  Rice,  "The  Mathematical  Analysis  of  Random  Noise",  Bell 
Sys.  Tech.  J.  23(1944),  282-332,  24(1945),  46-156. 

[84]  F.J.  Kingsbury,  "An  Experimental  Model  Study  of  Underwater 
Acoustic  Scattering  from  a  Wind-driven  Surface",  Ph.D. 
Dissertation,  Dept,  of  Engineering  and  Applied  Science, 

Yale  University,  New  Haven,  Conn.  1970. 


R-6 


[85]  A.K.  Fung  and  R.K.  Moore,  "The  Correlation  Function  in 

Kirchoff's  Method  of  Solution  of  Scattering  of  Waves  from 
Statistically  Rough  Surfaces",  J.  Geophys.  Res.,  Vol.  71, 
No.  12,  2939-2943,  1966. 


> 


I 


I 


I 


I 


» 


